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ABSTRACT 

This  paper  illustrates  the  use  of  graphical  analyses  by  pre¬ 
senting  the  solution  of  six  problems  in  the  fields  of  operational 
analyses,  mechanics,  and  engineering:  The  Jeep  Problem,  the  Range 
of  a  Fleet  of  Aircraft,  a  Beam  under  Combined  Compression  and  Trans¬ 
verse  Load,  the  Problem  of  Car  Replacements,  Determination  of  Bal¬ 
listic  Trajectory  Parameters,  and  the  Two-Magnetic-Reactor  Problem. 

The  purpose  of  this  paper  is  to  arouse  an  interest  in  a  meth¬ 
odology  which  is  further  enhanced  by  the  graphical  display  capa¬ 
bility  available  in  today’s  computers  with  all  its  potential  prob¬ 
lem  solving  flexibility.  The  examples  treated  in  this  paper  are  not 
the  stereotyped  problems  forming  the  usual  subject  of  textbooks  on 
graphical  methods  and,  in  that  sense,  should  prove  of  greater  in¬ 
terest  to  the  reader. 
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I.  INTRODUCTION 


Graphical  solutions  have  not  been  popular  in  the  United  States ; 
and  now,  with  the  widespread  use  of  computers,  it  may  seem  even  more 
a  remote  possibility  to  arouse  interest  in  a  skill  not  widely  employed. 
Or.  the  other  hand,  the  graphical  input  and  display  capability  and 
computer  processing  of  graphical  information  may  yet  foster  the  ac¬ 
ceptance  of  graphical  solutions  and  prove  to  be  a  tool  to  train 
graphical  visualization. 

This  training  is  as  important  as  the  training  in  any  other  math¬ 
ematical  symbolism,  without  which  the  shorthand  of  operations  and  re¬ 
lationships  would  remain  as  hieroglyphics  before  discovery  of  the 
Rosetta  stone.  A  quite  natural  question  would  be  to  inquire  as  to 
the  reasons  for  previous  lack  of  enthusiasm  for  graphical  solutions. 
That  may  have  been  due  to  many  causes  :  lack  of  training  in  the  use 
of  graphics;  the  hard-to-understand  disdain  of  engineering  students 
of  their  own  language — the  language  of  drawings;  the  preponderance 
of  analytically  minded  mathematicians  to  geometers,  etc.  The  princi¬ 
pal  cause,  however,  may  have  been  unfounded  fear  of  the  lack  of  ac¬ 
curacy  of  graphical  solutions.  I  say  "unfounded”  because  more  often 
than  not  the  accuracy  of  a  graphical  solution  is  amply  adequate,  even 
without  having  to  alibi  it  by  mentioning  the  underlying  assumptions, 
physical  constants,  and  other  factors  which  preclude  cur  precise 
knowledge  of  a  given  physical  phenomenon. 

The  virtue  of  a  graphical  approach  is  twofold:  a  diagram  often 
suggests  to  a  trained  mind  a  solution,  wut  in  all  cases  contains  a 
visual  interplay  of  variables  indicating  their  relative  importance. 

In  that  sense  alone,  I  feel  justified  in  quoting  Oliver  Heaviside — 
although  his  statement  was  unrelated  to  graphics --that  there  is  no 
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better  way  to  prove  a  fact  than  to  show  it  to  be  a  fact.  Here  the 
use  of  the  words  "to  show"  is  intended  to  convey  an  idea  of  a  picture, 
a  diagram;  but,  of  course,  to  a  trained  mind  a  shorthand  symbol  is 
just  as  clear. 

I  have  selected  from  my  experience  six  problems — one  as  recent 
as  this  paper,  same  going  back  to  my  early  work  in  aircraft  engineering. 
Each  of  the  six  example  problems  selected  was  chosen  to  demonstrate  an 
approach,  to  make  a  point,  and  to  illustrate  the  solutions  of  practi¬ 
cal  problems  encountered  in  practice. 

The  first  two  problems,  the  jeep  problem  and  the  range  of  a 
fleet  of  aircraft,  were  selected  because  their  analytical  solutions 
as  offered  in  mathematical  literature  did  not  suggest  that  equally 
accurate,  and  in  fact  rigorous,  solutions  were  possible  using  consid¬ 
erably  more  elementary  graphical  approach.  The  solution  to  both  of 
these  problems  is  approached  by  drawing  diagrams  showing  the  relation¬ 
ship  of  the  fuel  consumed  as  a  function  of  distance,  and  in  the  jeep 
problem,  the  direction  of  travel.  The  actual  solution  becomes  under¬ 
standable  by  proper  juxtaposition  of  lines  and  figures. 

The  third  problem,  that  of  a  beam  under  combined  compression 
and  transverse  load,  illustrates  the  reaction  of  a  mind  disposed  toward 
graphical  approach.  The  expression  of  the  bending  moment  for  a  beam 
under  combined  compression  and  transverse  load  had  been  derived  long 
before  I  found  it  necessary  to  use  it.  Also,  its  use  was  widespiaad 
in  those  days  because  airplanes  were  biplanes,  and  in  wing  construc¬ 
tion  one  used  routed  wooden  spars.  The  critical  design  points  were 
not  obvious;  to  verify  the  adequacy  of  design,  the  stress  analyst  used 
his  desk  calculator  to  arrive  at  the  bending  moments  and  shears  and 
physical  properties  of  the  spar  cross  section  to  verify  local  factor 
of  safety.  It  was  perhaps  natural  that  with  ny  preference  for  graphics, 
that  type  of  solution  suggested  itself  because  of  the  trigonometric 
nature  of  the  analytical  expression  of  the  bending  moment.  In  reality, 
of  course,  these  functions  represent  the  limiting  magnitude  of  appro¬ 
priate  infinite  series. 
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The  fourth  problem,  that  of  car  replacements,  once  more  offered 
an  illustration  of  how  by  drawing  a  diagram  showing  the  relationship 
of  a  number  of  cars,  their  ages,  arid  stipulated  change  of  their  mean 
age  after  a  specified  period,  the  solution  suggested  itself.  Of  course, 
the  problem  as  posed  and  its  solution  are  based  on  the  assumption  of 
linear  variation  of  car  ages  and  the  discarding  of  oldest  cars  to  ar¬ 
rive  at  a  new  mean  car  age.  The  finding  of  the  answer  requires  care¬ 
ful  examination  of  the  diagram  and  the  ability  to  recognize  unimportant 
inaccuracy  in  areas  being  balanced  to  find  the  location  of  the  solution 
line. 

The  fifth  problem,  the  determination  of  ballistic  trajectory 
parameters,  is  interesting  because  it  calls  for  the  application  of 
several  disciplines  and  because  it  demonstrates  that  often  the  fear 
of  a  graphical  solution's  not  being  accurate  enough  is  ill-founded. 

This  example  problem  shows  comparison  of  the  results  obtained  graph¬ 
ically,  analytically  by  the  use  of  a  desk  calculator,  and  finally 
those  arrived  at  from  a  computer  program.  It  is  significant  that 
even  the  apogee  velocity,  which  graphically  is  obtained  as  a  differ¬ 
ence  of  two  relatively  large  numbers,  differs  by  less  than  one  percent 
from  the  computer  answer.  This  is  amply  accurate,  considering  the 
fact  that  in  this  particular  case  the  assumption  of  a  trajectory  in 
vacuum  and  a  nonrotating  earth  were  considered  acceptable  for  the 
purpose  for  which  the  problem  had  to  be  solved.  The  additional  inter¬ 
est  of  the  solution  resides  in  the  use  of  logarithmic  coordinates  em¬ 
ployed  to  save  certain  arithmetic  calculations. 

The  twc-magnetic-reactor  problem,  the  sixth  problem,  is  an 
illustration  of  a  problem  whose  solution  was  entirely  intuitive.  I 
assume  that  unless  once  is  inclined  to  think  graphically  one  would 
have  difficulty  in  sensing  the  steps  which  led  one  to  the  eventual 
solution.  On  the  other  hand,  once  there  was  a  graphical  solution  its 
transformation  to  a  grapho-analytical  form,  a  simpler  and  more  elegant 
form,  became  clear. 
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Appendix  A  contains  my  extension  of  the  use  of  the  logarithmic 
polar  diagram  of  an  airplane  proposed  and  used  by  me  as  a  young  engi¬ 
neer.  It  is  reported  here  because  it  is  the  original  use  of  the  log¬ 
arithmic  scales  by  Mr.  Rith  that  started  me  on  their  use  and  interest 
in  graphical  solutions  in  some  of  my  analytical  work. 


II.  THE  JEEP  PROBLEM* 
N.J.  Fine,  Washington,  D.C. 


1.  Introduction.  The  problem  in  logistics  with  which  this  paper  deals  was 
proposed  to  the  author  by  Gail  Young  and  Ivan  Niven,  both  of  Purdue  Univer¬ 
sity,  in  the  latter  part  of  1945.  The  original  source  is  unknown  to  the  author.  At 
that  time  Niven  had  obtained  a  partial  solution  based  on  certain  assumptions. 
After  the  first  submission  of  this  paper  it  was  learned  that  L.  Alaoglu  had  also 
obtained  a  complete  solution.  He  mentioned  that  similar  problems  had  arisen 
in  air  transport  operations  in  the  China  theater.  It  has  also  been  suggested  that 
there  may  be  applications  to  Arctic  expeditions  and  interplanetary  travel.  This 
paper,  however,  will  confine  itself  to  the  lowly  jeep. 

Suppose  that  a  jeep  can  carry  a  maximum  load  of  n  gallons  of  gas  and  can 
travel  c  miles  per  gallon.  The  jeep  is  required  to  cross  a  desvrt  *  miles  wide.  Our 
problem  is  to  prescribe  a  method  for  making  the  journey  most  economically  and 
to  find  the  least  sufficient  amount  of  gas.  It  is  not  obvious  that  such  a  method 
exists,  and  it  would  be  more  exact  to  speak  of  the  greatest  lower  bound,  until  the 
existence  of  the  minimum  is  established. 

We  shall  assume  that  «  and  c  are  both  unity.  This  involves  no  loss  of  gen¬ 
erality;  it  is  equivalent  to  taking  as  our  unit  of  distance  nc,  the  number  of  miles 
that  the  jeep  can  travel  on  a  full  load. 

If  x  ^  1,  the  problem  is  trivial.  If  x  exceeds  1,  however,  gas  dumps  will  have 
to  be  established  at  various  points  along  the  way.  It  will  be  convenient  to  take 
the  path  of  the  jeep  along  the  positive  x-axis,  starting  at  x  and  ending  at  the 
origin.  The  gas  dumps  will  then  form  a  subdivision  <r  of  the  interval  (0,  x) : 

<r.  0  <  x,  <  x>  <  •  ••  <  *r  <  x, 

in  which  the  x<  denote  the  positions  of  the  dumps  (assumed  to  be  finite  in  num¬ 
ber).  If  s  is  any  non-negative  number  less  than  x,  the  subdivision  <r  induces  a 
subdivision  of  (0,  i)  by  deletion  of  all  the  stations  to  the  right  of  *.  There  will 
be  no  ambiguity  if  we  refer  to  this  induced  subdivision  by  the  same  symbol,  9. 
Other  subdivisions  will  be  denoted  by  9*.  and  so  forth.  If  all  the  stations 
(points  of  division)  of  9  are  contained  among  the  stations  of  cr\  we  shall  say  that 
9'  is  a  rtfi-.vnsni  of  9,  written  9'  <9. 

We  may  now  rephrase  our  problem.  Once  a  subdivision  is  fixed,  the  amount 
of  gas  required  is  still  a  function  of  the  method  of  establishing  and  employing  its 
stations.  W’c  shall  denote  by/(x,  9)  the  greatest  lower  bound  of  this  amount  for 
all  possible  methods,  and  by/(x)  the  greatest  lower  bound  of  /(x,  9)  for  all  pos¬ 
sible  subdivisions  9  Our  task  is  to  discover  the  form  of /(x). 

In  {2  we  introduce  the  standard  method  of  establishing  and  using  the  sta¬ 
tions  of  a  given  subdivision  9,  and  we  prove  that  this  method  is  at  least  as  eco¬ 
nomical  as  any  other.  This  enables  us  to  determi,'.e/(x,  r)  in  |J.  A  rather  sur¬ 
prising  application  of  the  standard  method  leads  to  the  result  (|4)  that  if 
9'  <9,  then  f(x,  9')  5/(x,  #).  In  |5  we  determine  criteria  for  non-improvement 
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by  refinement.  These  criteria  lead  us  to  the  construction  of  an  optimum  <r*  and 
to  the  explicit  representation  of/(x,  a*)  -/(*)  (§6).  In  §7  we  derive  a  simple  and 
accurate  asymptotic  formula  for  f(x).  The  last  section  is  devoted  to  a  few  re¬ 
marks,  including  a  comparison  of  the  exact  solution  with  the  result  obtained 
by  considering  the  stations  equally  spaced  (one  of  Niven’s  assumptions). 

2.  The  standard  method.  One  very  natural  method  of  employing  the 
stations  of  a  given  <r  is  to  build  up  the  stockpile  of  gasoline  at  x,  by  making  all  the 
trips  between  *  and  x,  before  gang  to  to  *,_i,  and  to  continue  in  this  way 
throughout  the  journey.  In  other  words,  once  we  go  beyond  any  station  x,  we 
never  return  to  the  preceding  one,  x,+i. 

Suppose  that,  by  some  other  method,  m  complete  round  trips  are  made 
starting  at  x,  followed  by  a  last,  one-way  trip  from  x  to  xr.  The  ith  one  of  the 
round  trips  consists  of  Ait  the  one-way  trip  from  x  toxP;  B<,  the  round  trip  start¬ 
ing  and  ending  at  x,;  and  C„  th?  return  trip  from  xr  to  x.  Let  g<  be  the  amount 
of  gas  in  the  jeep  at  the  start  .of-^,.  Since  2(x-xr)  is  the  amount  used  in  per¬ 
forming  trips  Ai  and  Ci,  the  amount  g,-2(x-xr)  plus  the  residue  of  the  pre¬ 
ceding  trips  is  sufficient  to  perform  £,.  If  we  replace  the  sequence  Ah  B\ ,  Cj, 
A s,  Bt,  Ci,  ,  A n,  1^*,  C*,  by  A ,,  Ci,  A j,  Ct,  *  *  *  ,  A n,  Cm,  iin-fi,  and 
deposit  at  x,  the  amount  g<-2(x-xf)  after  each  4<(t  =  l,  2,  •  •  •  ,  m),  and 
gm+i  —  (x— x,)  after  Am+u  we  shall  then  be  in  a  position  to  perform  all  the  £,  in 
exactly  the  same  order  as  before.  When  this  has  been  done,  the  final  configura¬ 
tion  will  not  have  been  altered  and  no* more  gas  will  have  been  used.  The  same 
reasoning  applies  to  all  the  trips  starting  at  x„  and  so,  by  induction,  the  stand¬ 
ard  method  is  established  as  being  at  least  as  economical  as  any  other.  Hence¬ 
forth  we  shall  assume  its  use. 

3.  Determination  of  /(x,  «r).  Now  we  suppose  that  there  is  given  a  sub¬ 
division 


0-:  0  -  Xo  <  X,  <  <  Xr  <  Xr+i  **  X, 

and  that/(xi_i,  <r)  has  already  been  determined.  Cleariy/(x«,  cr)  =0,  so  we  have 
the  initial  step  in  the  inductive  definition  of  /(x,,  a).  Let  ft,  be  the  number  of 
trips  to  be  made  from  x,  to  x,_j.  Obviously  ft,g  1.  No  gas  is  to  be  left  behind  at 
xi(  since  that  would  imply  waste.  Hence  the  difference  between  /(x,,  a)  and 
/(x,_j,  a)  must  be  accounted  for  by  the  amount  used  in  the  2k,  —  1  trips  between 
the  two  stations.  Writing  A,  »=x,—  x»_i.  we  have 

(1)  /(*..  *)  -  /(*.-»,  *)  -  (2ft,  -  1)A,  (i  -  1,  ■  •  •  .  r  +  1). 

We  must  now  determine  ft,.  The  maximum  amount  of  gas  that  can  be  trans¬ 
ported  on  each  of  the  first  ft,  —  1  trips  is  1  -  2A, ;  on  the  last,  1  -  A,,  since  there  is 
no  return.  The  total  must  not  be  less  than  the  amount  required  to  proceed  from 
x,_i.  Hence 

(2)  *,(1  -  2A,)  +  A,  £  /(*•— i,  *). 
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From  (1)  it  is  clear  that  the  number  of  trips  must  be  as  small  as  possible  so  that 

(3)  (*.  -  1)(1  -  2A<)  4-  A,  <  /(*«->,  <0. 

provided  that  k,>l.  It  is  easy  to  see  that  i»  =  l  for  x,£\,  and  that  f(x ,,  o)=x, 
in  this  case.  If  xt>\,  we  have  £<>1.  In  this  case,  (2)  and  (3)  determine  the 
integer  k,  uniquely,  and /(*<,  o)  is  then  obtained  from  (1).  We  remark  that  if  the 
equality 

tn(  1  -  2A()  4-  A,  =  a) 
holds  for  some  integer  m,  then  i k,  =  m. 

We  shall  now  derive  a  useful  relationship  between  k,  and  f(x ,,  a).  If  we  elim¬ 
inate  f(x,-h  o)  between  (1)  and  (2),  we  have,  for  all  /£  1, 

(4)  k,  £  /(*„  a). 

Similarly,  (1)  and  (3)  yield 

(5)  *«  -  1  <  /(*„  a), 

provided  that  fc(>l.  But  we  see  directly  that  (S)  is  also  valid  for  k,  =  l,  so  (4) 
and  (5)  hold  for  all  /SI.  If  we  define  {a}  as  the  least  integer  not  less  than  a, 
then  for  all  Z2;  1  we  may  write 

(6)  *,=  {/(*„,)}. 

Since  f(z,  a)  is  an  increasing  function  of  t, 

(2)  k,  S  k<+i,  (  ^  1. 

Summing  (1),  we  obtain 

f+i 

(8)  /(*,  o)  =*  £  (2  k,  —  1)A(. 

4-1 

4.  Refinements  of  subdivisions.  Let  o'  be  a  refinement  of  a.  The  quantity 
f(x,  o)  may  be  thought  of  as  the  result  obtained  by  applying  a  non-standard 
method  to  o',  namely,  passing  over  those  stations  of  a'  which  do  not  belong  to  a. 
It  follows  immediately  from  §2  that  /(x,  o')  is  not  greater  than  f(x,  o),  that  is, 

W  /(*.  ”')  £  /<*,  o)  if  o'  <  o. 

ff  (x (_j,  x,)  is  an  interval  of  a,  with  the  associated  parameter  k,,  and  if 
yi),  (y\,  Vj),  ■  •  •  ,  (»_i,  y, *>x,)  are  intervals  of  o',  with  parameters 
k',  k",  •  ■  ■  ,  k(p),  then 

(1°)  *<p>  -  l  f(*. .o') }  £  (/(*„  v)|  - 

Using  (7),  we  obtain 

(ID 

From  (8), 


J V  S  k"  $  ■  ■  ■  s  *<»)  s 


(12) 


/(*..  o')  -  /(*.-!.  o')  -  £  (2*<‘>  -  1)(*  -  Vi-l). 

»-l 

Equation  (1)  may  be  written  in  the  form 

(13)  /(*.,  <r)  -  /(*,_,,  a)  =  £  (2*.  -  l)(y(  -  *_,). 

1*1 

Subtracting  the  members  of  (12)  from  those  of  (13),  we  have 

(14)  /(*„  <r)  -  /(*<,  o')  =  /(*,_,,  o)  -  /(*,_,,  <r')  +  2  2  (*,  -  *«>)(*  -  %_,). 

\Vc  observe  that  all  the  differences  in  (14)  are  non-negative  F'om  this  we  deduce 
that  actual  improvement  by  refinement  takes  place  if  and  only  if  we  can  find  an 
interval  (x(_i,  x,)  and  an  integer  i  such  that  kt>kli).  By  (11),  this  is  equivalent 
to  k,>k’. 

5.  Properties  of  a*.  Our  problem  will  be  solved  if  we  can  find  a  subdivision 
a*  for  which 

(A)  /(*,  o')  g  /(*,  <r)  for  every  a. 

We  can  bring  to  bear  the  results  of  §4  by  proving  that  any  o*  which  satisfies 
(A)  also  satisfies  (B)  that  follows,  and  conversely. 

(E)  /(*,  o*)  =  /( x,  o')  for  every  o'  <  a*. 

Clearly,  (A)  and  (9)  imply  (B).  Conversely,  suppose  that  (B)  is  satisfied,  and 
let  a  be  any  subdivision  whatsoever.  We  choose  for  o'  the  common  refinem  ;nt  of 
<r  and  o*.  From  (B),  /(*,  o*)=f(x,  o')\  another  application  of  (9)  shows  that 
f(x,  o')  5/(x,  a).  Combining  these  we  obtain  (A). 

Using  the  criterion  established  at  the  end  of  §4,  we  find  that  (B)  is  equivalent 
to 

(C)  For  every  <  =  1,  2,  •  •  •,  r  +  l,  and  for  every  y  such  that  x(_i<y^xi, 

\f(y.  O)  =  {/(*..  o*)}. 

We  shall  now  show  that  (C)  is  equivalent  to  (D) : 

(I))  For  every  m  - 1,  2,  •  •  ■  ,  [/(*,  o’*)],  there  exists  an  integer  s  such  thatf(x„  o •) 
=  m. 

Suppose  first  that  (C)  fails  for  some  t  and  y.  Since  /(*,  e*)  is  strictly  increas¬ 
ing, 

(15)  /(*,-.,  o*)  <  f(y,  o*). 

By  the  definitions  of  k'  and  the  function  (  |, 

(16)  f(y,  #•)  £  V. 
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Since  (C)  fails,  k'  <k,,  that  is, 

(17)  k'$kt-  1. 

Finally,  since  k,-  j /(*<,  a*) }  cannot  exceed  /(*,,  <r*)  by  as  much  as  unity, 

(18)  k,  -  1  <  /(*„  a*)- 

Combining  (15),  (16),  (17),  and  (18),  we  see  that  the  integer  k'  lies  strictly  be¬ 
tween  f(xt-u  a*)  and  f(x,,  a*).  By  monotonidty,  k'  cannot  be  equal  to  f(xt,  a*) 
for  any  s,  so  (D)  fails. 

Conversely,  if  m  is  an  integer  satisfying 

(19)  /(*<- 1.  »*)<»<  /(*»,  »*), 

we  can  find  a  number  y  such  that 

(20)  f(y,  «r*)  -  m, 
and  (C)  is  violated.  To  prove  this,  set 

.  _  ,  i.O 


y  =  *<-i  + 


2m  —  1 


and  let  A=y— x,_t.  Clearly  A  is  positive,  and 

(22)  m(l  -  2A)  +  A  -  /(x«_,,  ,*)■ 

Referring  to  the  remark  in  §3,  we  see  that  m  is  the  parameter  associated  with 
the  interval  y),  and 

(23)  f(y,  «*)  -  /(xt-u  O  +  (2m  -  l)(y  -  x,_0  -  m, 

which  proves  (20).  This  completes  the  proof  that  (D)  is  equivalent  to  (C),  (B), 
and  (A). 

6.  Construction  of  a*.  It  is  now  almost  trivial  to  construct  a  <r*  satisfying 
(D)  and  therefore  (A).  Merely  choose  the  stations  x*  so  that  /(x,*,  <r*)**t. 
Clearly  this  can  be  done  for  /  =  1 .  Suppose  that  x{,  •  •  •  ,  x,_i*  have  been  found. 
We  must  determine  x*  by 

(24)  *,-{/(*,*.  ,*)} 

(25)  /(x,*,  *»)  -  /(x,*_,.  „•)  -  (2*,  -  l)(x,*  -  x,*_ ,)• 

The  left  member  of  (25)  equals  unity;  hence 

(26)  x,*  -  x,*_,  -  (21  -  1)"‘. 

Therefore 


1  1 

„.-,  +  T  +  T  + 


2t  -  1 
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It  is  easy  to  verify  that  (27)  leads  to  the  required  equations  4, =/(*,•,  = 

The  subivision  <r*  so  determined  evidently  satisfies  (D).  If  r  is  the  greatest  in¬ 
teger  for  which  xr*  does  not  exceed  x,  we  may  write 

(28)  f(x)  -  /(x,  O  -  r  +  (2r  +  l)(x  -  **), 

and 


(29)  0£x-  x*<  (2 r+  l)-‘. 

It  has  now  been  shown  that/(x),  which  represents  the  number  of  gallons 
needed  to  take  the  jeep  *  miles,  is  a  function  which  is  piecewise  linear  over  in¬ 
tervals  of  length  1,  1/3,  1/5,  1/7,  and  so  on,  the  slope  of  the  graph  over  the  nth 
interval  being  the  nth  odd  number  so  that  the  function  takes  consecutive  in¬ 
tegral  values  at  the  corner  points. 


7.  An  asymptotic  formula  tot  fix).  From  equation  (1)  it  is  possible  to  get  a 
rough  idea  about  the  order  of  magnitude  of fix).  We  have  approximately 

A/  -  (24  -  l)Ax, 

4  -  fix). 

Neglecting  the  -1  compared  with  4,  we  find 

A/ 

—  -  2 Ax 

/ 

log  /  «•  2x  +  Ci, 

/  »  C*e**. 

We  shall  not  attempt  to  make  these  heuristic  methods  precise,  but  shall  pro¬ 
ceed  directly  to  a  derivation  based  on  the  exact  solution  obtained  in  the  pre¬ 
ceding  section.  Let  us  define 


(30) 


S{r) 


1  1 

1  +— +—  + 
2  3 


1 

+  — 
r 


It  is  well  known  that  for  large  r, 

(31)  Sir)  -  log  r  +  C  +  ~  +  0  , 

where  C  is  Euler’s  constant,  .577  •  •  •  ,  and  the  C(l/r’)  denotes  an  error  term 
whose  absolute  value  does  not  exceed  a  certain  constant,  independent  of  r,  multi¬ 
plied  by  1/r’.  This  constant  may  vary  from  one  equation  to  the  next.  From  (27), 
(30),  and  (31),  we  have 

x.*  -  S(2r)  -  “5(r)  -  (log  (2r)  +  C  +  1  +  o(^j 
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Therefore, 


*r  =  log  (2y 7)  + 


Taking  exponentials, 


Squaring,  we  obtain 


“p  (** '  7 )  *  v,‘  (°(;;))  -  2V-(i  +o(-i)). 


~  exp  (2*f*  -  C)  -  r  +  0 


Thus,  writing  g(u)  =  J  exp  (2«  -  C), 

(35)  /(*,♦)  -  f  -  *(*•)  +  0  -  ,(*•)  +  0(e-«>. 

Now  for  any  m  satisfying 

/36)  **  <  «  *  *  <  x*+1, 

we  have,  by  equation  (28), 

(37)  /'(«)  «  2r  +  1. 

Also, 


(38)  «'(«)  "  2g(«)  >  2g(x*)  .  2r  +  0  (-1), 

(39)  «'(*)  <  2g(**+l)  »  2r  +  2  +  0 

(40)  /'(«)  -  *'(«)  =  0(1). 

Integrating  (40)  between  x  *  and  x,  we  find 

/(*)  ~  «(*)  -  (/(*,*)  -  !(**))  +  f’(/'(u)  -  g'(u))du 

m0(y)  +  o(*-  *?)  -o(- i) 

-  0(e~u). 

Therefore,  for  all  x, 

^  /(*)  ■  i  exp  (2*  -  C)  +  0(e-**). 

dMi'LbrbhTv'i  “l" "™r  "rm  i!  n°'  ”*>  th«n  the  pri„. 

P  .  but  it  actually  tenda  to  zero  exponentially.  Alzo,  the  approximate 
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equation  derived  by  such  high-handed  methods  at  the  beginning  of  this  section 
turns  out  \o  be  much  more  accurate  than  might  have  been  expected. 

8.  Remarks.  Suppose  we  decide  to  use  equal  subdivisions  of  the  interval 
(0,  x ),  N  in  number.  We  can  prove  that  if  N  is  chosen  to  vary  with  x  so  that 
x1/#  tends  to  zero,  then  the  increase  in  the  amount  of  gas  required  (over  the 
exact  solution)  is  less  than  50%  for  large  x.  To  counterbalance  this  increase, 
however,  the  number  of  stations  required  will  be  very  much  smaller.  For  exam¬ 
ple,  if  N  varies  as  the  cube  of  x,  it  will  also  vary  as  the  cube  of  log  r.  Hence,  if 
we  take  into  account  the  cost  (in  time,  energy,  material,  and  so  on)  of  setting 
up  the  stations,  it  might  very  well  happen  that  the  equal  subdivisions  would  be 
more  economical.  Of  course,  we  should  then  have  an  entirely  new  problem. 

For  a  fixed  *  and  sufficiently  large  N,  we  can  come  as  close  as  we  please  to 
the  minimum  /(x)  by  means  of  N  equal  subdivisions.  An  amusing  fact  here  is  that 
the  minimum  can  be  attained  if  and  only  if  x  is  rational. 

We  close  with  several  remarks  about  the  character  of  the  solution  a*  we  have 
obtained  here.  It  is  obvious  that  this  solution  is  not  unique,  since  every  refine¬ 
ment  of  cr*  is  also  a  solution.  It  can  be  shown  that  the  converse  is  also  true;  that 
is,  every  solution  is  a  refinement  of  a*.  Furthermore,  if  we  consider  the  class  of 
subdivisions  a  for  which  the  function  f(x,  a)  is  continuous,  we  can  prove  that  this 
class  is  identical  with  the  class  of  all  solutions.  It  would  be  of  interest  to  see 
whether  this  criterion  can  be  obtained  directly,  and  whether  the  minimum /(x) 
can  be  derived  from  it. 


12 


III.  GRAPHICAL  SOLUTION  TO  EXAMPLE  PROBLEM  NO.  1  (THE  JEEP  PROBLEM) 

Postulates  for  optimum  solution  are : 

1.  The  number  of  trips  between  two  consecutive  stations  is 
always  odd. 

2.  The  graph  of  the  fuel  consumed  during  the  entire  trip, 
viewed  from  the  terminus,  is  a  series  of  straight  lines 
having  a  continually  increasing  slope. 

3.  To  obtain  maximum  efficiency,  the  jeep  must  always  ini¬ 
tiate  each  trip  to  the  next  dump  station  fully  loaded. 

Consistent  with  Dr.  Fine’s  choice  for  n  and  c,  we  take  n  =  c  =  1. 
This  means  that  if  the  ordinate  axis  is  c  when  the  abscissa  axis  is  n, 
a  45  deg  line  gives  the  relationship  between  miles  traveled  and  fuel 
consumed. 

Figure  1  graphically  portrays  the  relationship  between  possible 
fuel  dump  locations  (c^),  fuel  consumed  in  traveling  between  each  pair 
of  proximate  points  (c  F  ,  F  F'),  and  fuel  available  after  a  single 
trip  between  each  such  pair.  The  word  "possible”  in  specifying  the 
fuel  dump  locations  and  the  figure  are  drawn  consonant  with  the  post¬ 
ulants  I  and  III  and  Dr.  Fine's  choice  for  n  and  c. 
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Examining  Fig.  1,  we  can  see  that  segments  FnF^  represent  the 
quantity  of  fuel  which  can  be  deposited  at  the  fuel  dump  and  still 
leave  enough  fuel  (F^F")  to  travel  back  for  the  next  fuel  load  pickup. 
Thus,  if  the  fuel  dump  is  located  1/3  of  the  range  from  start,  a 
round  trip  consumes  2/3  of  the  range,  resulting  in  the  deposit  at  the 
dump  being  1/3  of  the  range.  A  trip  to  the  dump  without  return  per¬ 
mits  deposit  of  2/3  of  the  range. 

Similarly,  if  the  distance  between  dumps  is  1/5  range,  the 
round  trip  permits  deposit  of  3/5  of  the  range,  while  the  last  single 
trip  makes  it  possible  to  deposit  4/5  of  the  range. 

The  efficiency  of  depositing  fuel,  i.e.,  the  ratio  of  the  de¬ 
posited  fuel  to  that  consumed  in  the  round  trip,  increases  as  the 
distance  between  fuel  dump  points  decreases. 


FIGURE  1. 


We  now  proceed  in  Fig.  2  to  construct  a  graphical  solution  of 
a  possible  relationship  between  fuel  requirements  and  locations  of 
fuel  dumps  for  ranges  exceeding  that  of  the  jeep.  Consistent  with 
postulate  III,  we  start  at  the  terminus  and  establish  fuel  dump  pre¬ 
vious  to  it  as  being  a  full  jeep  range  away.  The  triangle  OAF  estab¬ 
lishes  the  relationship  between  fuel  consumed  (the  ordinate)  and  dis¬ 
tance  traveled  (the  abscissa).  The  location  of  the  nearest  fuel  dump 
can  be  inferred  from  the  observation  that,  consistent  with  postulates 
I  and  III,  x  =  1/3  is  a  logical  contender  for  the  preferred  distance. 
T  '.e  figure  shows  the  three  trips  involved  in  the  transfer  and  estab¬ 
lishes  the  sufficiency  of  point  as  a  solution. 


Total 
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One  observes  now  that  the  next  fuel  dump  must  contain  sufficient 
amount  to  permit  the  deposit  cf  two  full  ranges*  at  fuel  dump  1/3. 

This  requires  preferably  increased  efficiency  of  fuel  transfer.  Thus, 
it  is  logical  to  assume  the  new  fuel  dump  location  somewhat  closer, 
which  suggests  the  point  1/5  as  the  choice.  In  Fig.  3  it  is  readily 
established  by  drawing  the  appropriate  45  deg  lines  representing  the 
trips  between  the  1/3  point  and  the  new  point  that  five  such  trips 
consume  one  jeep  range  while  affording  the  deposition  of  two  ranges 
at  point  1/3,  thus  producing  point  F^.  Using  similar  reasoning,  we 
construct  points  F^,  Fg,  etc.,  corresponding  to  progressively  de¬ 
creasing  distance  between  fuel  dumps. 

Thus,  we  have  arrived  at  a  solution  which  is  consonant  with  our 
three  postulates. 

Now  we  shall  endeavor  to  prove  that  this  solution  represents 
the  preferred  solution. 

The  1/3  range  distance  seems  to  be  so  logical  that  it  is  reason¬ 
able  to  question  the  introduction  of  the  1/5  and  1/7  fuel  dumps.  As 
a  test  we  will  eliminate  fuel  dump  1/5  and  replace  fuel  dump  1/7  with 
a  new  fuel  dump  located,  from  fuel  dump  1/3,  a  distance  of  another 
1/3  range.  Figure  4  shows  that  this  choice  of  a  new  fuel  Comp  penal¬ 
izes  by  one  full  range— from  F?  to  F' — not  even  considering  the  slight 
loss  of  the  distance. 


The  word  "range”  is  used  here  to  designate  the  amount  of  fuel 
needed  for  that  distance. 
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This  illustration  shows  that,  if  instead  cf 
starting  from  point  D  and  depositing  fuel 
at  C  and  proceeding  to  B,  one  started  from 
distant  C/3  from  B  and  wer.i  directly  to  B, 
this  would  entail  an  additional  full  range  of 
travel,  thus  requiring  initially  at  D  .  not  n  =  4 
but  n  -  5;  BD  =  '2/35  while  , 
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FIGURE  4. 
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Similar  tests  (Fig.  5)  will  readily  prove  that  fuel  dump  loca¬ 
tions  closer  than  those  originally  chosen  in  conformance  with  the 
stated  postulates  will  not  affect  the  fuel  requirements  but  will  in¬ 
crease  the  number  of  trips,  while  locations  greater  than  the  original 
ones  will  result  in  severe  penalty  of  the  required  fuel. 


This  illustration  proves  that  there  is  no 
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IV.  THE  RANGE  OF  A  FLEET  OF  AIRCRAFT*# 

J.  N.  FRANKLINf 

1.  Introduction.  The  problem  discussed  in  this  paper  is  to  determine  the 
range  of  a  fleet  of  n  aircraft  with  fuel  capacities  g,  gallons  and  fuel  effi¬ 
ciencies  ti  gallons  per  mile  (i  =  1,  •  •  •  ,  n).  It  is  assumed  that  the  aircraft 
may  share  fuel  in  flight  and  that  any  of  the  aircraft  may  be  abandoned  at 
any  stage.  The  range  is  defined  to  be  the  greatest  distance  which  can  be 
attained  in  this  way.  Initially  the  fleet  is  supposed  to  have  g  gallons  of 
fuel. 

A  theoretical  solution  is  obtained  by  the  method  which  Richard  Bell¬ 
man  [1]  calls  dynamic  programming.  Explicit  solutions  are  obtained  in  the 
case  of  two  aircraft  with  different  fuel  capacities  and  fuel  efficiencies  and 
in  the  case  of  any  number  of  aircraft  with  identical  fuel  capacities  and 
identical  fuel  efficiencies. 

The  problem  is  similar  to  the  so-called  jeep  problem.  The  jeep  problem 
was  solved  rigorously  by  N.  J.  Fine  [2],  A  solution  was  also  obtained  by  0. 
Helmer  [3,  4].  Fine  cited  an  unpublished  solution  by  L.  Alaoglu.  The  prob¬ 
lem  was  generalized  by  C.  G.  Phipps  [5].  Phipps  informally  developed  the 
special  result  which  is  deduced  in  §  4  of  this  paper. 

2.  A  recurrence  formula.  Let  Cm  be  any  subset  of  m  of  the  given  n  aircraft. 
Mathematically  Cm  may  be  represented  by  a  subset  of  m  of  the  first  n 
positive  integers.  Let  M(g,  Cm)  be  defined  as  the  range  of  the  fleet  of  m 
aircraft  C„  starting  with  g  gallons  of  fuel.  Then  the  required  range  of  the  n 
given  aircraft  is  M(g,  C„),  where  C,  =  1 1,  2,  ■  •  •  ,  n).  When  there  is  only 
one  aircraft, 

(1)  M(g,Ci)  =  min (g/r, ,  g,/r,), 

where  C|  consists  of  just  the  ith  aircraft. 

When  there  are  m  >  1  aircraft,  a  distance  x  is  flown  by  all  m  aircraft. 
Then  one  aircraft  is  abandoned,  leaving  a  subset  Cm-\  C  Cm .  It  is  un¬ 
necessary  to  consider  abandoning  more  than  one  aircraft  at  a  time.  For 
example,  the  effect  of  abandoning  two  aircraft  from  Cm  may  be  obtained 
by  abandoning  one  of  the  aircraft  from  Cm  and  then  immediately  aban¬ 
doning  the  second  aircraft  from  Cm-, .  After  the  distance  x  the  amount  of 
fuel  remaining  is 

(21  h  =  g  -  rL.tr„r,. 

The  greatest  distance  which  can  be  attained  by  the  remaining  aircraft 

*  Received  by  the  editors  Auguet  30, 1M6  end  in  revised  form  Mareh  17.  19fl0. 
t  California  Institute  of  Technology,  Paaadena,  California. 

#  Reproduced  by  permission  of  the  publisher. 
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Cm- 1  starting  with  h  gallons  of  fuel  is  M{h,  Cm_ i).  If  g,  Cm ,  and  x  are 
prescribed,  the  subset  Cm- 1  should  be  chosen  so  as  to  maximize  the  remain¬ 
ing  distance  M(h,  Cn- 1).  The  total  distance  traveled  wi’l  then  be 

(3)  X  +  m&Xcm-iClCmM (g  -  .r£cm  r, ,  Cm- 1). 


The  maximum  distance  M(g,  Cm)  is  obtained  by  maximizing  the  last 
expression  with  respect  to  x.  In  other  words,  if  g  is  £  the  total  capacity 
for  i  in  Cm  , 


(4) 


M(g,  Cm) 

=  max,S(,ic<l  „  [x  +  maxc„.lCcB  M(g 


X^,Cn  r,  ,  C  m-l  )  j- 


In  this  maximization  it  is  required  that  x  be  g!  £r,  i  8>nce  this  is  the 
greatest  distance  which  all  m  aircraft  esn  fly  with  g  gallons  of  fuel  before 
one  aircraft  is  abandoned.  If  g  is  given  in  excess  of  the  total  fuel  capacity 
Y^g,  of  the  aircraft  in  Cm  ,  then  some  fuel  must  be  thrown  away  and 

(5)  M(g,  Cm)  =  (0>£c„0,). 

The  recurrence  formulas  (4)  and  (5)  uniquely  determine  M(g,  Cm)  for 
all  subsets  Cm  with  m  =  2,  3,  •  •  •  ,  n.  It  is  easy  to  see  that  each  function 
M(g,  Cm)  is  polygonal  in  g,  i.e.,  continuous  and  piecewise  linear.  In  fact, 
(4)  may  be  rewritten,  by  the  identity  (2),  as 

(6)  M(g,  Cm)  =  maxkS(  [a(jf  -  h)  +  rnaxrm.lCr„  M(h,  C„_0|, 

where  a  -  l/£.cc»  r<  •  We  know  from  ( 1 )  that  every  function  M ( g ,  f, I 
is  polygonal.  Let  us  suppose  that  every  function  M(h,  Cm.x)  is  polygonal 
Then 

(7)  P{) 0  =  maxc„.,cc„  M(h,  CM_,) 


is  polygonal,  since  it  is  the  largest  of  a  finite  number  of  polygonal  fun< 
tions.  Now  (6)  takes  the  form 

(8)  M(g,  Cm)  =  ag  +  max*Sl)  \-ah  +  P{h)\. 

Since  —  ah  +  P(h )  is  polygonal,  its  maximum  value  for  h  £  g  is  a  poly- 
g(  nal  function  of  g,  say  P*(g),  and  therefore  M  =  ag  +  P*(g)  is  also 
polygonal  This  completes  an  inductive  proof  that  M{g,  Cm)  is  polygonal 
for  g  $  V( h  .  within  which  range  (8)  holds.  It  now  follows  from  (5)  that 
Mg,  Cm  is  polygonal  for  all  g.  Incidentally,  the  identity  (8)  shows  that 
Mig.  C,  i:  steadily  increasing  with  ra  e  So*  l/£>c<v  r,  when 
g  •:  Vj  ,r  'irger  values  of  g,  (5)  shows  that  M  is  constant. 
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3.  The  case  of  two  aircraft  ,  Lot 

Mlg)  =  M(y,  (\),  M i(g)  =  M(g,  C.)  (C,  =  t  =  1,  2). 

1  >>r  0  S  0i  +  </2,  (8)  takes  the  form 
(,0)  M(g)  =  ag  +  max*s„  [-ah  +  P(h)\, 

where 


<1^  P(h)  =  max, .1,2  A/, (A),  o  =  l/(r,  +  r2). 

By  the  identity  (1), 

(,2)  M,(h)  =  min  (A/r, ,  gx/n)  (i  =  1,2). 

Without  loss  of  generality  it  will  be  assumed  that  r,  g  r2 . 

Case  1.  Suppose  that  f/,/r,  ^  g2/r-t ,  In  this  ease 


HT)  l\h)  =  A/,(A)  (for  all  A). 

Then 

-ah  +  h/r,  ( A  g  y,) 

-aA  +  gCr,  (h  ^  0i). 

Si  nee  a  <  1/ri ,  it  follows  that 

(13)  /*•(,)  -  l-*" +  ^  (#S,,) 

[-agi  +  gjrx  (g  >  0,), 


(14)  -ah  +  l’(h)  = 


where  /J*(0)  =  max  [  —  ah  +  P(h)]  for  A 

[  o/ri 

(10)  (1/(0)  = 

la(ff  -  gi )  +  0,/r, 


g  0.  From  (8)  it  follows  that 
(g  ^  0i) 

(01  ^  0  s  0,  -F  0j). 


For  0  >  0i  +  02 ,  equation  (a)  gives 


(1 


M(g)  ■-  agt  +  g,/r ,  (0  >  gi  +  0i). 


From  the  definition  (2)  of  the  remaining  fuel  A  as  a  function  of  the  dis¬ 
tance  x  to  be  traveled  by  both  aircraft,  it  is  clear  that  the  optimal  procedure 
in  Case  I  is  to  use  just  aircraft  1  if  g  g  0i ,  or  if  0  >  0i  to  use  both  aircraft 
until  only  0i  gallons  of  fuel  remain  and  then  to  complete  the  trip  with  just 
aircraft  1, 

Case  2.  Suppose  that  0i/r,  <  gi/rt .  In  this  case 


(18) 


P(h)  - 


\  h/r:  (A  g  0i) 

0./r i  (01  S  h  g  0*) 

A/rj  (0*  £  h  g  0i) 

0?/rs  (A  g  0j) 
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where  g*  =  r^,/r,  is  the  abscissa  of  the  point  of  intersection  of  the  graph* 
of  Mi(h)  and  The  function  -ah  +  P(h)  is  a  polygonal  function 

with  peaks  at  h  -  Q\  and  h  =  gt .  There  are  two  subcases,  depending  upon 
whether  the  first  peak  is  higher  (A)  or  lower  (R)  than  the  second  peak. 

Case  BA.  Suppose  that  gx/rx  <  g,/r2  and  0,/r,"’  a  gt/r*.  Then 


(19) 


r*(g) 


-ag  +  g/r,  (g  g  gx) 
—agi  +  g,/r,  (g  £  g ,), 


and,  as  in  Case  1,  M{g)  has  the  form  ( 1(‘>),  ( 17).  The  optimal  procedure  in 
this  case  is  the  same  as  that  in  Case  1 . 

Cine  SB.  Suppose  that  <  gt/r,  and  <  y-./n.  Let  g'  =  gxr2/r? ; 
this  is  the  first  value  of  h  >  gx  at  which  —  ah  +  P(h)  =  ~agx  +  /'(</,). 
Then 


(20) 


P*(g) 


'  -ag  +  g/r\ 
-agi  +  gt/r, 
—  ag  +  ff/r5 
,  -ag-:  +  gt/r- 


(g  S  gt) 

(gi  £  g  S  g') 

(>/  £  g  S  g.,) 

(g  ^  gt). 


Therefore, 

(21)  Mg)  = 


g/r\ 

a(g  -  f/i)  +  gCn 
g/rt 

a(g  —  gt)  +  gt/r-t 


( g  s  9l) 

(ffi  ^  f/  £  s') 

(0'  £  p  £  0:) 

(?»  £  0  S  0i  +  5  ), 


and,  according  to  equation  (5), 

(22)  M(g)  =  ay ,  +  f /,  r.  (p  >  p,  f  </.  i 


The  optimal  procedure  is  as  follows.  If  y  g  ,  use  only  aircraft  1  If 
p,  g  p  g  p',  use  both  aircraft  until  only  pi  gallons  remain;  then  use  ju>t 
aircraft  l.  If  /  £  g  £  gt ,  use  only  aircraft  2.  If  g:  £  g  g  g,  +  gt ,  use 
both  aircraft  until  only  g .  gallons  remain;  then  use  just  aircraft  2.  If  <j 
>  ffi  +  Qt .  some  fuel  must  be  thrown  away,  and  the  trip  is  made  with 
g  =  (/,  +  g-,  as  descrilted  in  me  preceding  sentence. 

From  these  results  it  is  apparent  that  in  the  general  case  ol  n  aircraft 
the  optimal  policy  will  depend  in  a  complicated  way  upon  y  as  well 
upon  the  p,  and  r,  .  For  example,  the  value  of  y  may  determine  which  ol 
the  aircraft  finishes  the  trip. 
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4.  The  case  of  identical  aircraft.  Let 

* =  G,  r,  =  R  (i  =  1,  •••  ,  n), 

and  let  M m(g)  =  M(g,  Cm)  (m  =  !,•••,  n).  In  this  ease  the  recurrence 
formulas  (4),  (5)  take  the  form 

(24)  Mm(g)  =  max, s |.r  +  3/ m_,(g  -  mfta)l  ( g  £  mG), 

(2n)  Mjg)  =  J»/m(mG)  >  mG ). 

It  will  he  shown  that,  it  =  [ g/(!\,  the  greatest  integer  £  g/G, 

,w  "-(*>-S(,  +  3+-  +d  +  wrm 

(27)  Mm(g)  =  ^  ^1  +  ^  +  ’  ’  ’  +  ^  (0  >  mG). 

The  right-hand  side  of  (26)  is  defined  as  g/R  when  k  =  0.  In  the  optimal 
policy,  if  g  =  kG  £  mG,  the  trip  is  begun  with  k  aircraft.  If  kG  <  g 
<  (k  +  1  )G  5  mG,  the  trip  is  begun  with  k  -f  1  aircraft.  If  g  >  mG,  then 
g  —  mG  gallons  of  fuel  must  he  thrown  away,  and  the  trip  is  begun  with 
all  m  aircraft.  In  any  ease,  if  the  trip  is  begun  with  K  aircraft,  the  first 
aircraft  is  abandoned  whet  only  (K  —  1  )G  gallons  of  fuel  remain.  Then 
K  —  1  aircraft  are  flown  until  only  (K  —  2)G  gallons  remain,  and  so  on. 

This  result  can  he  e-ta  ilishcd  by  induction.  If  m  =  1,  formulas  (26) 
and  (27)  become 

(28)  -l/,(g)  =  g/R  (g  £  G), 

(2!D  Mx(g)  =  G/R  ( g>G ), 

which  is  correct  acc<  rd  ng  to  (1).  Assume  that  the  result  holds  for  m  -  1 
aircraft.  Then 

CIO)  ~  <l*)  =  7 - (0  <  h  <  (m  -  1)G) 

ah  (m  —  \)R 

at  all  points  h  at  '  Inch  the  polygonal  function  3/ _j( A )  has  a  derivative. 
Therefore,  x  +  .1/  ,tg  -  mRx)  is  a  steadily  decreasing  function  of  x  for 
0  <  j- <  g/mlt  if  /  S  (m  -  ])R.  Setting  x  ■*  0  in  (24)  gives  the  maximum 
value 

0<1>  Mm(g)  =  Mm-\(g)  * g  £  («  -  1K»>- 

Hut  .)/„  ,(g)i-  two  by  the  right-hand  side  of  (26)  for  ail  y  £  (m  -  1)G. 
Therefore,  (26 »  s  established  for  g  g  (m  -  MG. 

Next  suppo-i  that  (  »t  —  1  )G  <  g  £  mG.  Iheu  3/ ■.-■((/  —  mRx)  is 
constant  for  /  •  mRx  £  (m  -  1  )G;  for  larger  values  of  x  the  rate  of  in- 
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crease  is  g  —  mR/ (m  —  1  )R  <  —1.  Therefore,  the  maximum  (24)  is 
attained  when  g  —  mRx  =  (m  -  1  )G,  and 


(32) 


Mn( g)  =  °- - {~rl)G  +  -  1  )G) 


_  g  -  (m  -  1)G  ^ 
mR 


Gr(1+1+  ••  +^n)- 


This  establishes  the  result  (26)  for  (m  -  I)G  <  g  £  mG.  The  result  (27) 
for  g  >  mG  follows  from  (5). 

6.  An  asymptotic  formula  for  g.  The  solution  (26),  (27)  in  the  case  of 
identical  aircraft  is  similar  to  the  solution  of  the  jeep  problem,  although  the 
solutions  were  established  by  different  methods.  In  this  section  an  asymp¬ 
totic  formula  will  be  developed  for  the  amount  of  fuel  g  which  is  necessary 
in  order  to  transport  identical  aircraft  a  distance  .r.  Let  g  =  f(x).  It  will 
be  shown  that 


(33)  f(x)  =  -4(x)  +  0(exp  (  —  Rx/G)), 
where 

(34)  A(x)  =  (?(-i  +  exp(^-  C’)). 


In  these  identities  G  and  R  are  the  fuel  capacity  and  the  fuel  efficiency  of 
each  of  the  aircraft,  and  C  is  Euler’s  constant,  .577  -  .  This  result  is 

comparable  to  Fine’s  asymptotic  formula  for  the  solution  of  the  jeep  prob¬ 
lem  [2]. 

From  the  result  of  the  last  section  it  is  clear  that  x  is  the  range  of  n  +  1 
aircraft  with  initial  fuel  supply  g,  where  nG  <  g  £  (n  +  1)G.  Setting 
m  =  n  +  1  in  (26)  gives 


(35) 


x 


X*  + 


g  —  nG 
(n  +  1  )R' 


where 


(36) 


x* 


SO 

U 


+i  + 


+ 


0 


log  n  +  C  +  4  0 

£n 


This  well-known  asymptotic  formula  is  derived  in  [G,  p.  529).  From  (36)  i* 
follows  that 


(37)  „p(*|:-c)-ne,P(i  +  o(l)) 


»  +  j  +  0 
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This  relation  shows  that 


(38) 

(39) 


n  i? 


constant  -f-  exp 


/ftz* 
\  G 


'  =  0  (exp  (~Rx*/G), 


(40)  n  m  -i  +  exp  -  c)  +  0  (exp  (-fix VO). 
Since  nG  =  /(z*),  multiplication  by  G  gives 


(41)  fix*)  =  A(z*)  +  0(exp  (—Rx*/G)). 

In  order  to  justify  replacement  of  z*  by  z  in  the  identity  (41),  it  is  con- 
venient  first  to  show  that 


(42)  /'(z)  -  A’(x)  =  0(1)  (0  <  z  -  x*  <  G/(n  +  1)11). 

Differentiation  of  g  as  a  function  of  z  in  (35)  gives 

(43)  f(x)  =  (n  +  1  )R. 

But 

•4'(z*)  <  A’(x)  =  ft  ex pfe  -  c) 

(44)  \G  ) 

g  A'(x*  +  G/(n  +  DR)  g  A'(x*)  exp  (l/(n  +  1)). 
Since,  by  (39),  A'(x*)/(n  +  1)  *  0(1),  it  follows  that 

(45)  A'(x)  =  -4'(z*)  +  0(1). 

Subtraction  of  (45)  from  (43)  gives 

(46)  fix)  -  A'(x)  =  (»  +  DR  -  ft  exp  -  c)  +  0(1). 

The  required  relation  (42)  now  follows  from  (40).  Integration  of  (42) 
gives 

fix)  -  A(x)  =  fix*)  -  A(z*)  +  0(z  -  x*) 

(47) 

»  0(exp  (-  Rx*/G))  +  0(l/n)  -  0(exp(—  Rx*/G)). 


Since  z  —  x*  is  bounded,  x*  may  be  replaced  by  z  in  the  last  expression, 
and  this  gives  the  asymptotic  formula  (33). 
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V.  GRAPHICAL  SOLUTION  TO  EXAMPLE  NO.  2  (THE  RANGE 
OF  A  FLEET  OF  AIRCRAFT) 

I  am  indebted  to  my  colleague  Dr.  H.  Morris,  who  encouraged  me 
to  clarify  the  presentation  and  to  provide  a  treatment  of  the  case  of 
aircraft  of  unequal  speed.  Two  solutions  are  presented. 

A.  CASE  OF  AIRCRAFT  WITH  EQUAL  SPEEDS 

Initially,  the  problem  will  be  solved  using  the  two  tacit  assump¬ 
tions  contained  in  Franklin's  solution,  namely: 

1.  All  aircraft  are  flying  at  the  same  speed. 

2.  Time  lost  in  refueling  is  neglected. 
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B 

FIGURE  1. 


FIGURE  2. 


FIGURE  3. 


Max  range 


-J 


Figure  1  shows  the  relationship  between  fuel  consumed  and  dis¬ 
tance  flown;  for  example,  the  fuel  required  to  travel  distance  OB  is 
given  by  the  ordinate  OD.  The  relative  efficiencies  of  aircraft  are 
readily  indicated  by  the  slope  of  lines  AO- -the  smaller  the  slope,  the 
more  efficient  is  the  aircraft. 

Figure  2  shows  the  solution  for  the  case  of  two  aircraft  of 
equal  efficiency  and  equal  range.  In  this  case  it  is  immaterial  which 
of  the  two  aircraft  gives  up  its  fuel — the  maximum  range  remains  the 
same  in  either  case.  To  determine  the  range  and  location  for  refueling, 
for  convenience  of  graphical  solution  the  fuel-distance  relationship 
for  these  aircraft  is  drawn,  one  below  and  the  other  above  the  x-axis, 

OA  applying  to  aircraft  A  and  OB  to  aircraft  B. 

If  we  draw  line  A'C  parallel  to  line  OB,  the  ordinates  CC7  and 
C/C//  give  the  amounts  of  fuel  required  by  aircrafts  A  and  B,  respec¬ 
tively,  to  cover  the  distance  OC7.  It  becomes  evident  that  at  point 
C,  aircraft  B,  by  giving  up  all  its  fuel,  can  replenish  the  fuel  con¬ 
sumed  by  aircraft  A.  In  fact,  CC"  =  A/B,  and  CC7  -■  A  7B  -  C'C" ;  there¬ 
fore,  at  C7,  aircraft  B  has  just  enough  fuel  left  to  restore  the  fuel 
originally  carried  by  aircraft  A.  Thus,  aircraft  A  can  now  proceed 
to  a  new  destination  A"  such  that  C7A"  =  OA7. 

Next,  we  consider  the  case  of  two  aircraft  of  equal  efficiency 
but  unequal  range.  In  this  case  it  is  not  self-evident  which  of  the 
two  aircraft  should  give  up  its  fuel.  Examining  both  cases,  as  shown 
in  Fig.  3,  we  see  that  aircraft  B  should  refuel  aircraft  A.  The  graph¬ 
ical  construction  involved  is  the  same  as  before  but  is  repeated  for 
each  aircraft  to  establish  two  maximum  range  points:  point  B7  if  air¬ 
craft  B  is  used  to  refuel  aircraft  A  at  a  distance  OC7  from  start  of 
flight,  and  point  A7  if  aircraft  A  refuels  aircraft  B  at  a  distance 
0D#  from  the  departure  point. 


Figures  4,  5,  and  6  show  the  determination  of  maximum  ranges  for 


pairs  of  aircraft  of  differing  efficiencies  and  ranges. 
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FIGURE  6. 


We  will  assume  now  that  fleet  of  aircraft  examples  shown  in 
Figs.  4,  5,  and  6  represented  six  aircraft  flying  together.  Having 
established  the  preferred  order  for  refueling  for  each  pair,  we  are 
now  left  with  three  aircraft,  A^,  A^,  and  Bg,  the  subscript  designating 
the  figure  where  we  have  determined  the  new  range  of  the  respective  air¬ 
craft.  We  also  established  the  locations  where  the  refuelings  took 
place.  Figures  7  and  8  reproduce  these  locations  at  the  time  of  re¬ 
fueling,  also  individual  ranges  which  each  aircraft  could  cover  with¬ 
out  further  refueling.  For  clarity  we  redraw  Figs.  4,  5,  and  6,  com¬ 
bining  the  fuel-distance  relationship  of  the  three  aircraft  in  Fig.  7 
to  facilitate  graphical  solution  of  the  problem.  If  aircraft  B  were 
to  refuel  aircraft  A,  the  amount  of  available  fuel  is  given  by  line 
O'B^.  Proceeding  as  before,  we  find  point  A^  and  thus  the  range  ex¬ 
tension  of  aircraft  A^  to  point  A"  from  the  sequence  of  refueling  A'" 
with  aircraft  B  and  finally  with  aircraft  A,.. 

Figure  8  is  merely  a  check  to  see  whether  a  different  refueling 

sequence  would  improve  the  range.  In  the  illustration,  A4  refuels  A^, 

and  B,.  completes  the  refuleing  cylce  by  giving  up  its  fuel  to  Ac.  The 
b  j 

range,  although  close  to  that  found  in  Fig.  7,  is,  however,  smaller. 


B.  CASE  OF  AIRCRAFT  OF  DIFFERENT  SPEEDS  AND  EFFICIENCIES 

This  case  is  shown  in  Fig.  9.  In  addition  to  the  abscissae  of 
distance  we  show  to  the  left  of  the  origin  an  abscissa  of  time.  Thus, 
lines  OA^  and  OB^  give  the  time -fuel  relationship  for  aircraft  A  and  B, 
respectively.  Selecting  time  interval  t  from  start — assuming  in  this 
case  that  both  aircraft  depart  simultaneously — we  find  that  aircraft 
A  would  have  outdistanced  aircraft  B  by  G/A/.  For  aircraft  A  to  ar¬ 
rive  at  G7,  it  would  have  been  necessary  to  delay  its  departure 
by  a  segment  AC,  using  the  scale  of  time,  or  to  move  its  departure 
point  back  a  distance  AC,  using  the  scale  of  distances.  On  the  other 
hand,  if  we  wanted  to  have  the  two  aircraft  meet  at  a  point  A7,  it 
would  have  been  necessary  for  aircraft  B  to  be  stationed  distance  G7A7 
ahead  of  aircraft  A  or  to  have  delayed  the  departure  of  aircraft  A  by 
a  time  interval  represented  by  a  segment  G7A7,  using  the  scale  of  times. 

Because  of  linear  relationship  of  the  variables  involved  it  can 
be  seen  that  lines  OC  and  OD  represent  the  locus  of  points  through 
which  pass  lines  OA  and  OB  appropriate  to  selected  distances  of  ren¬ 
dezvous  of  the  two  aircraft.  Thus,  line  CC7  would  describe  a  condition 
in  which  aircraft  A  would  arrive  at  point  G7  when  aircraft  B  has  con¬ 
sumed  G7G  amount  of  fuel.  This  suggests  that  if  we  drew  lines  B'F 
parallel  to  line  OB,  and  B7I  parallel  to  line  OA,  we  would  establish 
points  H  and  E,  giving  the  earliest  time  for  complete  refueling  of 
aircraft  A  by  aircraft  B,  and  aircraft  B  by  aircraft  A,  respectively. 

By  inspection  it  is  evident  that  the  latter  choice  results  in  better 
range.  The  new  range  is  at  point  B'". 

It  is  believed  that  this  method  can  be  readily  extended  to  the 
case  of  several  aircraft.  Although  it  requires  a  trial-and -error 
solution  because  it  is  not  evident  in  what  combinations  refueling 
should  proceed,  the  simplicity  of  the  solution  results  in  relatively 
small  labor. 
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VI.  GRAPHICAL  SOLUTION  TO  EXAMPLE  NO.  3  (A  BEAM  UNDER 
COMBINED  COMPRESSION  AND  TRANSVERSE  LOAD) 

A.  PRECISE  BENDING  MOMENT 

The  general  expression  of  the  bending  moment  can  be  written  as 

M  =  Cx  Sin  (x/j)  +  C2  Cos  (x/j)  +  wj2  (1) 

where  x  is  the  distance  from  the  left  support  and  w  is  the  uniformly 
distributed  transverse  load.  In  the  expression  of  j  =  ,  E  is  the 

modulus  of  elasticity,  I  the  effective  moment  of  inertia  and  P  the 
axial  compression.  C^  and  C2  are  the  constants  of  integration  which 
can  be  determined  from  the  conditions  at  the  two  supports;  C-^  = 

[D2  "  D1  Cos  CL/ji )  5  C2  =  D1;  D1  and  D2  are  introduced  for 

brevity  and  designate  as  follows:  D1  =  M1  -  wj2;  D2  =  “  WJ2>  L  is 

the  length  of  the  bay,  the  bending  moment  at  the  left  support,  M2 
is  the  moment  at  the  right  support.  See  Fig,  1. 

Attention  is  invited  to  the  mathematical  significance  of  the 
trigonometric  functions  appearing  in  the  solution  of  differential 
equations  involved  in  the  problem  of  combined  compression  and  bending. 

These  trigonometric  functions  must  be  thought  of  actually  as  rep¬ 
resenting  the  infinite  series  whose  limiting  magnitude  can  be  considered 
as  the  sides  of  a  right  triangle  usually  employed  to  define  the  Sine  x 
and  Cosine  x,  when  one  angle  of  a  right  triangle  is  x  radians  in  mag¬ 
nitude.  This  property  makes  possible  the  graphical  construction  des¬ 
cribed  in  this  article,  but  the  reader  should  not  lose  sight  of  the 
true  significance  of  the  functions. 
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FIGURE  1. 


While  the  present  solution  can  be  extended  to  include  the  case 

of  transverse  loading  other  than  uniform  the  present  article  deals 

with  a  uniformly  distributed  transverse  loading  only.  Examining 

equation  1,  we  observe  that  by  drawing  two  circles  with  radii  equal 

o 

respectively  to  and  C2  and  distant  wj  and  by  considering  the  signs 
oi  the  products  Sin  (x/j)  and  Cos  (x/j),  it  is  possible  easily 
to  obtain  their  magnitude  for  any  value  of  x  from  zero  to  L.  It  will 
be  found  in  most  practical  cases  that  C,  and  C_  will  have  a  sign  oppo- 

n  -1.  £ 

site  to  that  of  wj''1  and  we  will  assume  that  they  are  negative.  The 
construction,  however,  is  simple  and  once  it  is  understood  it  can  be 
easily  changed  to  take  care  of  the  actual  signs  of  and  .  Noticing 
that  Sin  (x/j)  changes  from  zero  to  Sin  (L/j)  while  Cos  (x/j)  from  1 
to  Cos  (L/j)  it  is  obvious  that  the  angle  57.3  deg  (L/j)  must  be  drawn 
in  such  a  way  (if  we  want  to  measure  the  bending  moment,  along  y  axis) 
that  for  the  value  of  x  =  o  the  products  C^  Sin  (x/j)  and  0,^  Cos  (x/j) 
become  zero  and  C,,,  respectively.  This  can  be  obtained  by  measuring 
the  angle  57.3  deg  (L/j)  on  the  circle  C^  from  the  horizontal,  while 
on  the  circle  C^  it  must  be  measured  from  a  vertical  diameter. 


Draw  the  beam  1-2  and  on  the  continuation  of  the  line  choose  a 

center  from  which  the  circle  Cn  must  be  described.  Measure  downward 

2 

along  the  continuation  of  the  vertical  diameter  the  distance  wj  and 
from  that  point  as  a  center  describe  the  circle  C^.  Lay  out  the  angle 
equal  to  57.3  deg  (L/j)  from  a  horizontal  diameter  and  construct  this 
angle  also  on  the  circle  C2,  in  which  case  it  must  be  started  from  the 
vertical  diameter.  (The  construction  is  illustrated  in  Fig.  2  and 
made  clockwise.)  Dividing  the  length  of  the  bay  and  the  arcs  1-2  of 
the  circles  C-^  and  into  any  number  of  equal  parts,  we  construct  the 
bending  moment  diagram  by  simply  projecting  the  points  horizontally 
from  the  circles  on  the  vertical  lines  passing  through  the  corresponding 
points  on  the  beam.  In  Fig.  2  the  construction  is  illustrated  for  the 
point  b.  It  can  be  easily  proved  that  the  ordinates  enclosed  by  the  lines 
lb,  2^  and  1^  b2  2^  represent  the  bending  moments  at  the  corresponding 
points  of  the  beam  in  the  scale  in  which  were  drawn  the  circles  C,  and 
C2<  In  fact  b1  b2  =  Ob-^  +  0b2  =  Sin  (x/j)  -  wj2  +  C2  Cos  (x/j). 


FIGURE  2.  Construction  of  Precise  Sending  Moment  Diagram 


In  order  to  obtain  the  circles  C^  and  C it  is,  however,  not 

necessary  to  calculate  the  magnitude  of  C^  and  C2,  since  we  can  make 

the  following  observations:  From  Fig.  3,  we  see  that  since  the  dis- 

0 

tance  0  02  is  equal  to  wj  ,  if  we  lay  out  downward  from  0^  the  seg¬ 

ment  01 j  such  that  it  is  equal  to  we  obtain  immediately  C2> 

C2  =  D1  =  M1  -  wj2  =  -  (wj2  -  Mx). 

2  1 

To  obtain  the  =  lay  out  uPwarc^  ^2  eQual  t0  M2> 

obtaining  0^  0^  which  is  equal  to  -D2-  To  obtain  ,  after  having 
constructed  the  angles  57.3  deg  (L/j)  we  observe  from  Fig.  3  that 
taking  O2  2^  and  adding  it  downward  to  along  vertical  diameter  we 
obtain  point  2^  such  that  the  intersection  of  a  horizontal  through 
that  point  with  radius  0^  2^  gives  us  C.. 

02  2V  =  D1  Cos  (L/j) 

°1  =  d2  +  D1  Cos  (L/J} 

°1  21  =  °1  2y/Sin(L/j)  =  -  [D2  -  D1Cos(L/j)]/Sin(L/j). 

The  attention  is  invited  to  the  fact  that  angle  57.3  deg  (L/j)  must 
be  drawn  accurately  since  otherwise  it  may  result  in  a  large  error  in 
C,.  As  a  check,  the  distance  2^,  equal  to  M0 ,  may  be  used  before 
proceeding  with  further  construction. 

B.  PRECISE  SHEAR 

By  differentiating  the  expression  of  the  bending  moments  (equa¬ 
tion  .1)  we  obtain  the  expression  for  shear  as  follows: 

d.M  _  1  ,  _  _  x  ,  c.  x. 

dx  j  1  1  2  y 

It  can  be  seen  that  if  we  draw  the  beam,  divided  into  the  same  number 
of  parts  as  before,  along  the  continuation  of  the  vertical  diameter  and 

project  from  the  circles  C^  and  C2  the  points  1  . b  ........ 2  on 

the  horizontal  lines  drawn  through  the  corresponding  points  on  the  beam, 


the  horizontal  distances  between  the  two  curves  thus  obtained  will 
give  us  the  magnitude  of  shear  in  the  scale  of  M/j  where  M  is  the 
scale  of  moments  (Fig.  4). 
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FIGURE  3.  Circles  and  C ^ 

2 

FIGURE  4.  Construction  of  Precise 
Shear  Diagram 

C.  MAXIMUM  BENDING  MOMENT 

The  location  of  the  maximum  bending  moment  is  determined  from 
the  shear  diagram  at  the  point  of  zero  shear. 

D.  DEFLECTIONS 

To  obtain  the  deflection,  we  observe  the  fact  that  the  expression 
for  precise  bending  moment  can  be  written 

M  =  Mo  -  by  (2) 

Where  is  the  primary  bending  moment  due  to  the  transverse  load  alone. 
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Relationship  of  construction  points  to  structure  of  the  beom 


Scale  of  moments:  1"”1 2,000  in,  lb 
Scale  o*  deflections:  1"  -  P/M 
-6000/12,000  *  1 /2" 


Precise  bending  moment  diagram:  1"  “  20,000  in.  lb 
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Graphical  Solution  of  the  Beam  Analyzed  in  NACA  Technical  Note  383 
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FIGURE  5. 


From  equation  2  we  can  express  y  =  (MQ  -  M)/P. 

It  is  evident,  therefore,  that  in  order  to  obtain  the  deflection,  we 
can  simply  superimpose  the  diagram  of  precise  and  primary  moments  and 
by  dividing  the  difference  between  the  two  moments  by  axial  load  P,  to 
determine  the  magnitude  of  the  deflection  at  any  point.  If,  however, 
ve  should  replot  the  precise  bending  moment  diagiam  on  a  horizontal 
base  by  using  a  new  scale  of  moments  M  =  P  (where  M  is  the  new  scale 
cf  moments  and  P  is  the  compression  load)  and  construct  the  primary 
moment  to  the  same  scale,  we  will  obtain  directly  the  deflections  to 
their  actual  magnitude.  The  replotting  of  the  precise  bending  moment 
can  be  easily  accomplished  by  using  a  proportional  divider  or  by  a 
method  of  similar  triangles.  The  primary  bending  moment  for  a  uni¬ 
formly  distributed  load  is  simply  constructed  as  a  parabola  and  the 
example  (Fig.  5)  gives  the  constructional  lir  involved. 

E„  EXAMPLE  SHOWING  APPLICATION  OF  THE  METHOD 

To  illustrate  the  practical  application  of  the  method  described 
in  this  article,  a  complete  determination  of  precise  bending  moment, 
shear  and  deflection  diagrams  4 s  given  here  for  the  metal  truss  beam 
whose  properties  and  effective  moment  of  inertia  were  determined  by 
Mr.  Andrew  E.  Swichard  in  NACA  T.N.  393  (pages  24  to  31  inclusive). 

Data : 

L  =  180  inches:  P  =  -  6000  pounds;  I  eff  =6.76  inches4;  E  =  10,000,000 
pounds/square  :h;  M  =  -  40,000  inch  pounds;  M  0;  w  =  -  12  pounds/ 
inches;  j  =  97.98;  j2  =  9640;  w j x  =  -  115,200;  L/j  1.838  =  105  deg 
191;  wL2/8  =  48,600  inch  pounds. 


VII.  GRAPHICAL  SOLUTION  TO  EXAMPLE  NO.  4  (THE  PROBLEM 

OF  CAR  REPLACEMENTS) 

A  given  railroad  owns  N  number  of  cars,  purchased  over  a  period 

of  T  years.  It  is  assumed  that  the  ages  of  the  cars  can  be  expressed 

T 

as  a  linear  function  of  time,  thus  giving  ^  as  the  average  age  of  the 
cars.  It  is  proposed  to  reduce  the  average  age  to  a  stipulated  figure 
by  purchasing  new  cars  during  the  next  a T  period.  Assuming  that  the 
replacements  are  being  accomplished  by  yearly  purchases  and  that  car 
ages  can  be  considered  a  linear  function  of  time,  determine  the  number 
of  cars  to  be  purchased  and  replaced  in  order  to  maintain  the  original 
number  of  cars  but  lower  the  average  age  of  all  the  cars  at  the  end  of 
AT  period  to  the  desired  figure. 

In  Fig.  1,  at  start,  CD  is  the  average  age  of  the  N  cars.  The 
average  age  is  reduced  to  EF  after  the  required  number  of  cars  was 
retired  and  replaced  with  new  cars  during  the  interval  of  time 
AT  =  US  -  00  . 
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Assume  -  Linear  variation  of  cars'  ages. 
Time  and  cars'  ages  are  giver,  by  abscissa 
while  number  of  cars  by  the  ordinate 

CD  is  the  average  age  of  cld  cars 

EF  <  CD  average  age  of  cars  after 
replacement  of  oldest  cars  with  new 
cars  purchased  during  time  period  AT 

EF  is  stipulated  initially 


■j— Average- 
age  after  time 
I  interval  AT 


FIGURE  1 


Figure  2  shows  the  old  and  the  new  cars  at  the  end  of  the  AT 
replacement  period. 

Line  01  gives  the  ages  of  M  new  cars  purchased  during  the  period 
of  aT  =  years  (an  equal  number  of  old  cars  was  retired);  at  the 
end  of  the  00^  period  the  ages  of  the  remaining  old  cars  are  given  by 
line  AC.  By  averaging  the  ages  of  all  cars  at  the  end  of  the  00^  per¬ 
iod,  and  drawing  a  straight  line  BEL  passing  through  point  E,  such 
that  EF  is  the  specified  new  mean  age  of  all  the  cars,  we  define  an 
area  LBO^  gi/ing  the  total  age  of  new  and  remaining  old  cars. 

This  implies  by  inspection  that  area  AGHB  must  be  equal  to  area 
LHIO  or,  which  is  the  same,  that  area  APLB  is  equal  to  GPOI. 


FIGURE  2. 


One  can  resort  to  a  graphical  determination  of  M  (Fig.  3)  by 
proceeding  as  follows : 

Draw  line  RK  equidistant  from  lines  AO  and  NO^.  Erect  an  ordi¬ 
nate  at  L  intersecting  lines  PA  at  S.  Continue  PA  until  it  intersects 
KR  at  K'.  Mark  off  an  ordinate  K"R  equal  to  the  segment  K'T'  on  the 
line  KR.  By  inspection  we  state  that  with  sufficient  accuracy,  area 


Number  of  cars 


LSK'T'  is  equal  to  the  area  LSK"R.  If  at  midpoint  of  line  3K"  we 
draw  line  GI  parallel  to  the  abscissa  axis,  and  a  line  joining  I  with 
0,  the  figure  PGIO  gives  the  required  area,  while  its  ordinate  10  de 
fines  the  number  of  cars,  answering  the  problem. 

Figure  4  gives  an  example  of  the  application  of  this  method. 
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FIGURE  4. 


VIII.  GRAPHICAL  SOLUTION  TO  EXAMPLE  NO.  5  (DETERMINATION 
OF  BALLISTIC  TRAJECTORY  PARAMETERS) 

For  much  of  the  conceptual  and  analytical  work  dealing  with 
ballistic  missiles,  computer  precision  is  not  required. 

Given  two  of  the  three  basic  parameters  (burnout  angle,  burnout 
speed,  and  range  or  range  angle),  the  method  described  in  this  paper 
can  be  used  to  determine  vacuum  trajectory,  velocity  at  any  point, 
and  time  of  flight,  using  only  a  slide  rule,  compass,  and  ruler.  Figure 
1,  taken  from  "Free  Flight  of  a  Ballistic  Missile,"  by  Albert  D.  Whelan 
(American  Rocket  Society  Journal,  December  1959),  gives  the  relation¬ 
ship  of  burnout  speed,  burnout  angle  (from  vertical),  and  range  angle. 
The  method  assumes  a  nonrotating  earth  and  applies  to  the  portion  of 
the  unpowered  trajectory  outside  the  atmosphere.  For  many  problems 
the  data  are  adequate  for  a  complete  trajectory  from  ground  launch  to 
reentry. 


FIGURE  1.  Relationship  of  Burnout  Speed,  Burnout  Angle,  and  the  Range  Angle 


A.  BALLISTIC  TRAJECTORY  CONSTRUCTION 

Graphical  construction  of  the  ballistic  trajectory  is  made  on 
the  basis  of  the  following  laws  of  mechanics  and  properties  of  conics: 

1.  Ballistic  trajectory  is  an  ellipse  whose  one  locus  is  the 
center  of  the  earth. 

2.  Included  angles  between  a  tangent  to  the  ellipse  and  lines 
joining  the  point  of  tangency  with  each  focus  are  equal. 

2.  Sum  of  the  length  of  the  two  lines  joining  any  point  on  the 
ellipse  and  its  two  foci  is  constant. 

4.  Given  its  major  and  minor  axis,  an  ellipse  can  be  constructd 
by  the  use  of  the  auxiliary  circle. 

Figure  2  shows  the  construction  of  ballistic  trajectory  ellipse. 


FIGURE  2.  Construction  of  Ballistic  Trajectory  Ellipse 


B.  DETERMINATION  OF  SPEEDS  AND  TIMES  ALONG  THE  TRAJECTORY 


To  determine  missile  speeds  and  times  along  the  trajectory  we 
make  use  of  the  energy  relationship  which  gives  the  interplay  between 
kinetic  and  potential  energies : 
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where  VQ  is  the  burnout  speed;  V  is  missile  speed  at  some  point  N;  gQ 
is  gravitational  constant  at  sea  level;  is  a  segment  giving  the 
height  of  the  missile  from  the  center  of  the  earth. 

After  integration  we  obtain 
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Here  H  is  the  altitude  of  the  missile,  =  Ro+h,  and  V  is  the  speed 
at  that  point. 


From  Equation  2  it  can  be  seen  that  previously  obtained  trajec¬ 
tory  can  be  transformed  into  the  curve  of  potential  energy  by  correcting 

F-o 

each  altitude  by  a  factor  of  as  shown  in  Fig.  3. 


FIGURE  3.  Determination  of  Kinetic  and  Potential  Energy 

It  then  becomes  apparent  that  if  NN^,  the  segment  representing  poten¬ 
tial  energy,  is  subtracted  from  the  segment  representing  the  initial 
kinetic  energy,  LL^  =  NN^,  the  difference  thus  obtained,  repre¬ 

sents,  to  a  proper  scale,  the  vel'city  of  the  missile  at  that  point 
on  the  trajectory.  Assuming  that  we  have  obtained  such  a  segment, 
then  if  speeds  are  taken  in  feet  per  second,  accelerations  ii.  feet 
per  second  squared,  and  distances  in  nautical  miles,  the  speeds  would 
be  obtained  by  measuring  the  length  of  the  segment  using  the  scale  of 
distances  and  then  taking  a  square  root  of  the  product ; 


(Segment  length)  x  2  gQ  x  6080,  so  that  missile  velocity  at 
any  point  is 


V  =  625  ,/Segment  length  =  625 

If,  now,  we  were  to  divide  the  trajectory  into  a  number  of  equal 
segments  and  determine  at  midpoint  of  each  such  segment  the  average 
velocity  V,  then  the  ratio  of  the  segment  length  divided  by  the  aver¬ 
age  velocity  would  give  us  the  time  of  flight  along  that  segment  of 
trajectory. 

C.  USE  OF  LOGARITHMIC  SCALES  FOR  DETERMINATION  OF  SPEEDS  AND  TIMES 

The  graphical  construction  and  calculation  of  speeds  and  times 
can  be  more  conveniently  accomplished  by  the  use  of  logarithmic  scales. 

In  Fig.  4,  if  we  make  lgVQ  =  OL,  we  will  have  a  convenient  loga¬ 
rithmic  scale  (2)  of  velocities  and  At. 

Having  previously  divided  the  trajectory  into  a  number  of  equal 
lengt-h  segments  L-l,  1-2,  2-3,  etc.,  one  measures  lengths  N^N^  with 
the  linear  scale  of  lengths  1  which  in  Fig.  4  we  plotted  to  the  left 
of  0.  Next,  using  a  logarithmic  scale,  these  magnitudes  of  these 
length^  are  plotted  on  cor'  spending  rays  as  an  extension  of  lg  625 
which  was  marked  as  an  arc  of  a  circle  with  0  as  its  center.  Because 
the  magnitude  required  is  ^  lg  N"Nj',  the  logarithmic  scale  (3)  of  half 
the  magnitude  of  that  used  for  lg  625  must  be  employed.  A  smooth  curve 
LL^  plotted  through  these  points  gives  at  any  point  on  the  trajectory 
the  speed  at  that  point. 


Ti  rf'  TMFQ 
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Using  logarithmic  scale  (2)  ,  we  draw  an  arc  of  a  circle  with 
center  at  Fj  and  radius  equal  to  lg  (400x6080)  because  the  length  of 
the  equal  segments  into  which  we  divided  the  trajectory  is  400  nmi. 
Thus , 
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FIGURE  4.  Determination  of  Speeds,  Times  and  Plot  of  Complete  Trajectory  Data 


Ig  At  =  lg  ^|  =  ig  —&=  =  ^  -  Off  625  +  1  Ig  HJHJ), 
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expressed  in  terms  of  logarithmic  scale  (2) ,  gives  logarithm  of  At 
at  the  appropriate  points  1,  2,  3,  etc.  These  lengths  terminate  on 
the  arc  L  ^ne  Lm  and  the  speed  curve  LL^  as  shown  in  Fig.  4,  By  sum¬ 
mation  of  the  AtTs,  one  obtains  the  elapsed  times  from  launch  to  any 
point  on  the  trajectory. 

E.  EXAMPLE 

Figure  4  was  drawn  for  a  case  in  which  the  given  initial  condi¬ 
tions  are 


VQ  =  24,000  ft /sec 

Yq  =  41  deg  (from  vertical)  . 

From  Fig.  1,  one  finds  the  angle  9  =67.8  deg,  equivalent  ro  a 
range  of  4070  nmi.  Figure  4  gives  the  appropriate  scales  necessary 
to  construct  and  to  determine  missile  altitude,  velocity,  and  time  of 
flight  at  any  point  on  the  trajectory  and  summarizes  the  data  graphi¬ 
cally.  The  apogee  is  found  to  be  1370  nmi,  the  apogee  speed  10,160 
ft /sec,  and  the  total  time  of  flight  46.9  min.  The  toi~al  time  of 
flight  was  obtained  as  the  sum  of  At's  for  each  segment  as  shown  in 
the  following  tabulation. 
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Seconds 


No. 

Seqment 

At 

SAt 

1 

L-l 

118 

118 

2 

1-2 

130 

248 

3 

2-3 

146 

394 

4 

3-4 

163 

557 

5 

4-5 

184 

741 

6 

5-6 

205 

946 

7 

6-7 

221 

1167 

8 

7-8 

244 

1411 

9 

8-7 

244 

1655 

10 

7-6 

221 

1876 

11 

6-5 

205 

2081 

12 

5-4 

184 

2265 

13 

4-3 

163 

2428 

14 

3-2 

146 

2574 

15 

2-1 

130 

2704 

16 

1-L 

118 

2822 

I  am  indebted  to  Dr.  R.  Finke,  who  calculated  principal  traiec 
tory  parameters,  using  both  a  desk  calculator  and  slide  rule  as  well 
as  the  computer. 

The  following  tabulation  gives  the  three  sets  of  figures  and, 
thus,  a  relative  precision  of  the  graphical  and  analytical  methods: 


Graphical  Analytical  Computer 


Range,  deg 

68 

67.74 

67.74 

Range,  nmi 

4,070 

4,071.7 

4,072.66 

Apogee  altitude, 

nmi 

1,870 

1,858.1 

1,857.38 

Apogee  velocity, 

ft /sec 

10,160 

10,225.3 

10,228.81 

Flight  time,  min 

47 

46.78 

46.78 
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IX. 


GRAPHICAL  SOLUTION  TO  EXAMPLE  NO.  6  (THE  TWO-MAGNETIC 
REACTOR  PROBLEM) 


Given:  Two  magnetic  circuits  which  are  so  constructed  that 
they  share  some  volume  in  which  their  respective  magnetic  fields  and 
induction  (also  flux)  are  at  right  angles  (Fig.  1).  Each  core  may 
be  of  any  reasonable  material. 

Problem:  For  a  given  .4ttNI  in  the  control  magnetic  circuit, 
what  is  the  flux  through  the  power  circuit  for  various  ,4nNI  of  power. 


FIGURE  1. 
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A.  METHOD  A:  GRAPHICAL  SOLUTION 


Legend 

I  Current  in  amperes. 

N  Number  of  turns. 

I  Mean  length  of  the  magnetic  path. 

B  Flux  density  in  gauss. 

H  Magnetizing  force  in  oersteds  (gilberts/cm) 

F  Magnetomotive  force  in  gilberts. 

C,c  Properties  and  characteristics  of  the  control  reactor 
circuit. 

P,p  Properties  and  characteristics  of  the  power  reactor 
and  circuits. 

Z,z  Properties  and  characteristics  of  the  common  zone  of 
C  and  P  (Fig.  1). 

ZC&ZP  Subscript  denotes  characteristics  due  to  or  induced  by  the 
the  control  and  power  reactor,  respectively. 

An  examination  of  Fig.  1  suggests  the  following  relationships: 


Bzc  = 

Bc 

(1) 

BZP  = 

bp 

(2) 

HZC  +  HZP  ” 

Hz 

(3) 

2  2 

4 +  4  = 

Bz 

(4) 

In  addition,  since  it  is  possible  to  orient  the  common  zone 
material  in  such  a  way  as  to  ensure  identical  magnetization  properties 
at  right  angles  along  the  direction  of  core  C  and  core  P,  one  can 
write 


®ZC  _ 
^ZC 


(5) 
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The  magnetomotive  forces  in  the  control  and  the  power  circuit  are 
Fc  =  4H(t)dt  =  H  +  Hzc-tz  =  .4ttNcIc...  (6) 

Fp  =  »H(t)dt  =  Hp'tp  +  Hzpiz  -  .4irNpIp...  (7) 

The  relationships  B  vs  H  are  given  in  the  form  of  usual  mag¬ 
netization  curves,  indicated  schematically  in  Figs.  2,  3,  and  4. 


P 

FIGURE  3. 


1+- 14-0*25  * 

FIGURE  4. 

Also  known  are  the  number  of  tu^s  in  the  control  and  power 
reactor. 
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The  solution  of  the  problem  requires  the  determination  of  the 
magnitudes  of  H^,  Hp,  Hzc,  and  Hzp  for  any  selected  values  of  Fc 
and.  Fp. 

We  will  proceed  to  develop  a  graphical  solution  based  on  the 
following  observations; 


Assume  that  for  a  specific  value  of  and  a  specific  value  of 


Fp  we  obtained  the  corresponding  magnitudes  of  H^,  Hp,  Hzc,  and  Hzp 
2  2  2 

Then,  since  +  Hzp  -  Hz  we  can  spot  the  proper  points  on  the 
three  magnetization  curves  (as  shown  in  Figs.  5,  6,  and  7): 


t2 


We  redraw  curves  of  vs  and  Bp  vs  Hp  by  modifying  their 
abscissae  by  a  factor  —  and  j—  ,  respectively. 

This  is  shown  schematically  in  Figs.  8  and  9. 


B 


Note  that  the  difference  of  abscissae  in  these  curves,  namely 
FC 

q-Hoq  and 

Since  we  assumed  that  points  C  and  P,  and  hence  points  C1  and  P1, 
represent  the  solution  of  our  specific  problem  then  angles  B  in  Figs. 

10  and  11  are  equal  because,  in  that  case, 

which  fulfills  the  relation 

ships  (5)  and  (1)  and  (2). 


\ 

I 


A  check  of  the  correctness  of  our  solution  would  be  to  construct 
right  angle  triangle  shown  in  Fig*  12  below  and  compare  the  magnitude 
of  thus  obtained  with  the  Hz  of  Fig.  7. 


FIGURE  12. 


The  examination  of  Figs  10  and  11  suggests  a  graphical  solution 

of  the  problem.  We  plot  vs  Hc  j—  and  Bp  vs  Hp  having 

Z  Z 

selected  suitable  scales  for  the  coordinates. 

F  F 
C  P 

On  this  diagram  we  draw  vertical  lines  spaced  j—  and  —  from 

Z  h 

the  origin  of  coordinates  as  it  was  initially  shown  in  Figs.  10  and  11. 
On  curves  C  and  P  we  mark  several  pairs  of  points  corresponding  to  the 
number  of  equal  0^  angles  randomly  picked  by  us.  For  each  point  we 
obtain  a  pair  of  magnitudes  of  and  H^,  and  a  corresponding  Bp  and 
H—  among  which  accidently  there  may  be  a  point  representing  the  solu- 
tion  of  the  problem. 


6f 


Having  determined  we  obtain  Bc  and  Bp  as  well  as  and 
as  it  will  be  shown  later. 


The  solution  must  be  repeated  for  various  values  of  Fn  and  Fp 
and  cross-plotted  to  show  the  interdependence  of  the  selected  variables. 


B.  APPLICATION  OF  THE  METHOD 

^P 

Curves  B^  vs  H^  and  Bp  vs  Hp  £—  are  superimposed  on  the 
Z  Z 

same  graph  using  cross-section  paper. 

This  graph  is  used  as  a  template,  the  actual  construction  being 
performed  for  each  pair  of  Fc  and  Fp  on  a  separate  piece  of  tracing 
paper  which  thus  will  contain  all  details  of  construction  permitting 
an  easy  check  of  individual  solutions. 

For  one  group  of  solutions  either  F„  or  Fp  is  kept  constant. 

The  curve  corresponding  to  the  variable  factor  is  traced  on  the 
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I 

individual  vellum  and  placed  so  as  to  superimpose  vertical  lines  F  y 

z 

described  on  page  64.  With  this  arrangement,  a  gi’oup  of  radial 
lines  emanating  from  point  0^  provides  a  group  of  points  lc  and  lp, 

2c  and  2p,  3c  and  3p,  etc.  corresponding  to  points  P"1"  and  C’*'  of  Fig.  13. 


To  ensure  uniformity  of  all  solutions  we  will  agree  to  retain 
the  orientation  of  coordinates  of  curve  C,  for  example,  and  rotate 
90  deg  the  coordinates  of  curve  P. 

Having  marked  off  a  series  of  points  .l,  2,  3,  4,  etc.,  forming 
the  same  angle  with  the  coordinate  axes,  we  proceed  to  construct 
curve  8^  vs  H^. 

For  clarity  the  construction  is  shown  in  Fig.  16  for  point  2. 
Construction  of  vs  curve  illustrated  for  point  2.  Bc  and  Hzc 
ordinates  remain  oriented  as  drawn  vrhile  BD  and  H_„  orc'  .nate  rotated 

11  r  Zr 

90  deg.  Points  "Zn  define  2^  vs  curve. 
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After  repeating  this  construction  for  all  the  points  (normally 
4  to  5  points  suffice)  curve  vs  is  drawn  (Fig.  17)  and  super- 

vs  Hj,  curve  which  thus  determines  point 
Z  establishing  the  values  of  and  Bz  forming  the  solution  of 
the  problem.  The  determinat-*  on  of  point  Z  on  B„  vs  H„  magnetization 

u  a 

curve  is  shown  in  Fig.  18. 


imposed  on  the  selected  B^ 
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B 


Final  determination  of  all  the  unknowns  of  the  original  problem 
is  given  in  Fig.  19.  Having  obtained  point  Z  we  project  it  on  o'B' 
and  O'H^  axes  indicated  here  as  points  b£  and  H*.  The  magnitudes 
O'B"  and  O'H"  are  measured  off  on  O'B"  and  O'H"  curves  and  B_,  BD, 

L  u  L  r 

H^p  determined  directly  as  shown. 

The  practical  application  of  the  method  permits  a  number  of 
shortcuts : 

The  curves  need  not  be  redrawn  for  separate  constructions  and 
many  of  the  construction  lines  need  not  be  drawn  because  of  the 
reference  offered  by  the  graph  paper  on  which  are  drawn  the  magnetization 
curves. 
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C.  METHOD  B:  GRAPHO-ANALYTICAL  SOLUTION 

In  describing  the  alternate  method,  reference  will  be  made  to 
the  previous  explanation  of  the  graphical  solution  of  the  two-mag¬ 
netic-reactor  problem. 

In  this  case  the  magnetization  curves  are  drawn  on  log  log  paper 
having  introduced  on  this  paper  two  additional  scales  equal  to  one 
half  of  the  original  scales.  The  new  scales  thus  permit  direct 
reading  of  the  magnitude  of  the  square  of  any  ordinate  plotted  using 
the  original  scale  (Fig.  20). 

Instead  of  obtaining  graphically  the  curves  of  H„n  vs  E-  and 

ZtL  C 

H^p  vs  Bp  as  it  was  shown  in  Figs.  10  and  11,  a  table  shown  below 
is  prepared. 


.  .  i  •  1  •  I  •  I  < 


Curves  of  vs  Bp  and  H^p  vs  Bp  are  plotted  on  a  sheet  of  vellum, 
using  vrhe  previous  log  log  plot  as  a  template  as  shown  in  Fig.  21. 
Note  now  that  randomly  drawn  set  of  45  deg  lines  accomplishes  the 
same  purpose  as  the  radial  lines  of  Fig.  15.  In  fact,  because  the 
abscissae  and  c.dinate  scales  are  equal,  we  can  write 

log  Bq  -  log  Bp  =  log  Hzc  -  log  Hzp  or 
BC  _ 

b;  "  hzp  ••• 

which  fulfills  condition  (5). 
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FIGURE  21. 


Having  thus  obtained  a  series  of  lc  and  lp,  2c  and  2p,  3c  and  3p, 

etc.,  points  we  read  the  square  of  the  corresponding  ordinates,  enter 

them  into  Table  2,  add  the  respective  magnitudes,  and  plot  the  curve 
1  1 

of  B„  vs  H,,  whose  graphical  determination  was  given  in  Fig.  17.  In 
this  instance  the  operation  of  extracting  the  square  root  in  expressions 

BZ  =  VBzc  +  BZP  anc3  HZ  =VHZC  +  ^ZP  ta^en  care  using 
the  proper  scale. 


TABLE  2 


R2  *  R/2 

BC  Bzc 

B2  -  b'2 
BP  BZP 

R/2  ,  «/2 

BZC  *  BZP 

CM  O 
*3^ 

H/2 

HZP 

u#2  „/2 

HZC  HZP 

1 

1 

1 

i 

1 

» 

1 

i 

1 

i 

! 

Figure  22  shows  the  plot  of  B^  vs  H^  whose  intersect  with  the 
selected  Bz  vs  Hz  magnetization  curve  determines  the  point  Z. 

The  determination  of  points  C  and  P  and  corresponding  magnitudes 
of  B_j  B_,  H_,  H„-,,  and  H  is  given  in  the  illustration  and  is  self- 

C.  F  L  Li\*  L\Z 

explanatory. 


Example : 

Assumed  magnitudes 

In  f  p 

*r—  —  r—  =  5;  j*—  =  16.25 ;  ~  10 

Lz  ^Z  z 

Curve  C  is  Cobalt  Iron  (50*  Co,  50%  Fe) 

Curve  P  is  Medium  Silicon  Steel 
Curve  Z  is  Permalloy. 

Table  3  is  a  completed  version  of  Table  1,  and  Table  4  is  a  com¬ 
pleted  version  of  Table  2.  Figure  22  which  was  used  as  an  illustration 
of  the  method  also  represents  a  graphical  solution  to  the  values  con¬ 
tained  in  Tables  3  and  4, 
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FIGURE  22 


TABLE  3 


Bc 

Hc 

Hcr 

Hzc 

BP 

HP 

tp 

%  rz 

HZP 

1.77 

1.0 

5.0 

11.25 

1.3 

0.5 

2.5 

7.5 

4. 

1.5 

7.5 

8.75 

3.5 

0.75 

3.75 

6.5 

7.6 

2.0 

10.0 

6.25 

6.0 

1.0 

5.0 

5.0 

10.2 

2.5 

12.5 

3.75 

7.3 

1.25 

6.25 

3.75 

12.1 

3.0 

15.0 

]  .25 

8.4 

1.5 

7.5 

2.5 

13.3 

3.5 

17.5 

9.1 

1.75 

8.75 

1.25 

14.3 

4.0 

20.0 

9.8 

2.0 

10.0 

0 

15.5 

5.0 

25.0 

10.6 

2.5 

12.5 

16.5 

6.0 

30.0 

11.3 

3.0 

15.0 

17.2 

7.0 

35.0 

11.8 

3.5 

17.5 

17.7 

8.0 

40.0 

12.2 

4.0 

20.0 

18.2 

9.0 

45.0 

12.7 

5.0 

25.0 

18.7 

10.0 

50.0 

13.0 

6.0 

30.0 

20.0 

15.0 

13.3 

7.0 

35.0 

20.8 

20.0 

13.5 

8.0 

40.0 

21.4 

25.0 

13.7 

9.0 

4^.0 

21.8 

30.0 

13.8 

10.0 

50.0 

14.3 

15.0 

14.6 

20.0 

14.9 

30.0 

TABLE  4 


HZC 

8.7 

6.75 

10,5 

5.6 

11.9 

4.68 

12.8 

3.82 

14.1 

2.76 

H'2 

Hzc 

H'2 

HZP 

45.6 

30.0 

31.4 

19.6 

21.4 

13.6 

14.6 

9.0 

mm 

X.  CONCLUSIONS 

The  range  of  problems  of  the  six  examples  presented  here  is 
broad  enough  to  acquaint  a  reader  not  versed  in  the  subject  with  the 
variety  of  approaches  which  could  be  used  to  solve  a  given  problem. 

In  all  cases,  however,  graphical  problem  solving  has  one  common 
ingredient — to  prepare  at  the  outset  a  diagramatic  presentation  of 
the  relationship  of  variables  in  several  possible  forms.  In  fact, 
very  often  the  approach  of  classical  Euclidean  geometry  will  offer  a 
hint  of  required  construction.  By  the  classical  approach  I  mean  the 
assumption  that  the  problem  has  been  solved,  which  thus  permits  the 
drawing  of  the  solution  figure.  Inspection  of  this  figure,  with  its 
juxtaposed  variables,  will  often  disclose  their  principal  dependence 
or  relationship  pointing  to  a  solution.  A  good  illustration  of  this 
is  given  by  Figs.  2  and  3  in  the  problem  of  Car  Replacements  (pp.  47,48). 

The  example  problems  contained  in  this  paper  could  also  form 
instructive  exercises  in  training  one  in  the  use  of  graphical  computer 
display  and  techniques  of  interaction  between  the  computer  and  the 
operator. 

If  this  paper  will  whet  the  graphical  appetite  of  some  of  its 
readers  it  will  amply  justify  its  publication. 


APPENDIX  A.  A  NEW  APPLICATION  OF  THE  LOGARITHMIC  POLAR  DIAGRAM 

The  original  adaption  of  the  logarithmic  polar  diagram  is 
credited  to  Mr.  Rith  of  the  Eiffel  Laboratory  and  may  be  found  in 
Eiffel's  works  "La  Resistance  de  I’air  et  l’aviation"  and  "Nouvelles 
Recherches  sur  la  Resistance  de  I’air  et  l’aviation. "  Disadvantages 
of  the  original  Rith  and  later  methods  based  on  the  use  of  the  log¬ 
arithmic  polar  diagram  are  in  the  approximate  nature  of  the  estimated 
performance  obtained  and  the  consequent  inaccuracies  and  discrepancies 
as  compared  with  other  available  methods. 

THEORY  OF  THE  LOGARITHMIC  POLAR  DIAGRAM 

Rith  points  out  that  if  the  characteristic  curves  of  the  air- 
plai  e  are  plotted  as  a  polar  in  logarithmic  coordinates,  the  plot 
represents  a  functional  dependence  of  all  the  main  factors  entering 
into  the  two  principal  equations  of  the  flight  characteristics  of  an 
airplane : 

W  =  LV2  and 
HP  =  DV3; 

or  expressed  logarithmically: 

log  W  =  log  L  +  2  log  V  and 

log  HP  =  log  D  +  3  log  V  ; 

the ref ore : 

log  L  =  log  W  -  2  log  V  and 
log  D  =  log  HP  -  3  log  V. 

From  the  last  two  equations  it  is  apparent  that  we  can  consider 
log  L  as  the  sum  of  log  W  and  -2  log  V,  while  log  D  can  be  considered 
as  the  sum  of  log  HP  ar.d  -3  log  V.  The  values  of  log  W  and  log  HP 
can  be  plotted  on  the  direction  of  the  ordinate  and  abscissae  directly, 
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while  log  V  can  be  plotted  on  a  line  inclined  in  such  a  way  that  the 
values  projected  on  the  direction  of  the  ordinate  will  give  -2  log  V 
and  the  same  point  projected  on  the  axis  of  abscissae  will  give  a 
magnitude  of  -3  log  V.  This  new  axis,  called  the  speed  axis,  has  the 
slope  of  2/3  and  the  graduations  have  a  modulus  of  -  3.605 

times  that  of  the  modulus  used  on  the  other  axes. 

The  diagram  (Fig.  1)  thus  constructed  gives  the  necessary 
functional  dependence  of  the  factors  entering  into  the  two  equations. 


FIGURE  1 .  FIGURE  2.  Determination  of  Maximum 

Speed  at  Sea  Level 

In  order  to  take  into  account  the  effect  of  the  density  of  the 
air  and  so  to  obtain  the  characteristics  of  flight  at  altitude,  we 
observe  that  in  the  two  equations  mentioned  above  we  must  incorporate 
the  value  of  relative  density.  The  new  equations  are  written  as: 

W  =  n  L  V2  and 

HP  =  n  D  V3 

where  n  is  the  ratio  of  the  density  at  altitude  to  that  at  sea  level. 
The  neccessary  transformation  can  be  obtained  by  assuming  the  origin 
of  coordinates  moved  along  a  45-deg  line  so  that  the  magnitudes 
of  L  and  D  assume  a  new  value  equal  to  n  L  and  n  D.  Consequently, 
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in  order  to  obtain  the  relationship  of  the  factors  involved  at  altitude, 
it  is  sufficient  to  incorporate  a  new  reference  axis  (the  altitude 
axis)  which  has  a  slope  of  1/1,  and  to  plot  the  logarithmic  values  of 
relative  density  with  a  modulus  equal  to  a/I+T  =  1.41  of  the  modulus 
used  for  L  and  D  axes. 

Knowing  the  theory  of  the  logarithmic  polar  chart,  we  can  turn 
our  attention  to  the  practical  application  cf  it  to  performance 
estimate . 

It  is  assumed  that  the  characteristics  of  the  airplane  (that  is, 
the  lift  and  drag  at  one  mph  on  a  full-size  machine)  are  determined 
either  through  calculations  or  from  a  wind  tunnel  test.  Knowing 
the  values  of  L  and  D  and  having  a  sheet  of  coordinates  available, 
we  plot  the  logarithmic  polar  curve. 

DETERMINATION  OF  MAXIMUM  SPEED  AT  SEA  LEVEL 

On  the  right  hand  side  of  the  diagram  (Fig.  2),  a  vertical 
logarithmic  scale  is  marked  scale  |i.  This  scale  facilitates  the 
determination  of  the  available  thrust  horsepower  when  the  propeller 
efficiency  is  known.  Draw  a  horizontal  line  through  the  known 
propeller  efficiency  (line  1)  until  it  intersects  a  45-deg  line 
passing  through  the  point  of  BHP  maximum  (line  2).  Draw  a  vertical 
line  (3)  through  the  point  thus  obtained  until  it  intersects  a 
horizontal  line  (4)  passing  through  the  ordinate  indicating  the 
gross  flying  weight  of  the  airplane.  Point  A  is  the  thrust  horse¬ 
power  available  at  the  maximum  speed.  Drawing  through  point  B 
line  (5)  parallel  to  the  speed  axis,  the  segment  BC  obtained  between 
the  intersection  of  this  line  with  the  polar  (point  C)  and  point 
B  gives  the  magnitude  of  the  maximum  speed  of  the  airplane.  To 
find  its  value  lay  out  from  the  origin  of  coordinates  an  equal  segment 
on  the  speed  axis  and  the  reading  at  the  point  D  will  give  the  desired 
magnitude . 
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DETERMINATION  OF  STALLING  SPEED 


Draw  a  horizontal  line  (1)  tangent  to  the  uppermost  point  of 
the  polar  (Fig.  3).  Draw  line  (2)  through  given  gross  flying  weight. 
Line  (2)  is  parallel  to  the  speed  axis.  Segment  AW  (where  A  is  the 
point  of  intersection  of  the  two  lines)  gives  the  stalling  speed 
of  the  airplane,  which  must  be  measured  on  the  speed  axis. 


SPEED  FOR  FLIGHT  AT  MINIMUM  POWER 

Draw  line  (1)  tangent  to  the  polar  and  parallel  to  the  V  axis 
(Fig.  4).  Draw  line  (2)  through  W  corresponding  to  the  given  gross 
flying  weight.  The  intersection  of  the  two  lines  defines  point  A 
in  such  a  way  that  AB  is  the  speed  corresponding  to  minimum  power 
and  AW  is  the  minimum  power. 


FIGURE  3,  Determination  of  Stalling 
Speed  at  Sea  level 

MAXIMUM  L/D 
It 

mically 


FIGURE  4,  Determination  of  Speed 
for  Flight  at  Minimum 
Power 


vertical  scale  graduated  logarith- 


can  be  easily  proved  that  a 
with  the  same  modulus  as  used  for  the  L  and  D  axes,  and  located 
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arbitrarily  on  the  drawing  will  give  us  the  ratio  of  L/D  at  the  point 
of  intersection  of  this  scale  with  a  line  passing  through  the  cnosen 
point  and  inclined  at  an  angle  of  45  deg.  In  graduating  this  scale, 
we  must  determine  at  least  one  reference  point  and  graduate  the  scale 
in  the  same  sense  as  the  L  axis.  To  obtain  the  value  of  maximum 
L/D,  simply  draw  a  line  inclined  at  45  deg  and  tangent  to  the  polar 
(Fig.  5).  Point  A  gives  the  maximum  L/D  of  the  design. 

The  brief  explanation  of  the  theory  and  a  few  illustrations  of 
the  possible  uses  of  the  logarithmic  polar  diagram  given  above 
will  suffice  to  explain  the  practical  use  of  the  new  grapho-analytical 
method  which  consists  of  incorporating  the  curves  of  horsepower 
available  and  required  drawn  in  logarithmic  coordinates.  This  involves 
an  incorporation  of  a  few  more  reference  curves  as  explained  and 
illustrated  in  the  following  paragraphs. 

HORSEPOWER  REQUIRED  AT  SEA  LEVEL 

Knowing  V  max.,  calculate  values  of  0.9,  0.8,  0.7,  0.6,  0.5, 

0.4,  0.3  of  V  max.,  and  mark  them  on  the  speed  axis  as  shown  in 
Fig.  3.  This  is  not  absolutely  necessary  for  the  purpose  of  obtaining 
the  horsepower  required,  but  later  we  will  see  that  it  will  aid  us 
in  obtaining  the  horsepower  available.  The  drawing  of  the  curve  of 
horsepower  required  will  be  illustrated  by  following  the  procedure 
for  one  point;  for  instance,  point  V  -  0.5  V  max.  (Fig.  6). 

Draw  line  (1)  through  W  and  parallel  to  the  speed  axis.  Project 
point  0.5  V  max.  on  it.  This  gives  point  A.  Draw  a  horizontal  line 
through  this  point.  Through  the  point  B  draw  a  line  parallel  to 
speed  axis  until  it  intersects  the  horizontal  line  passing  through 
W.  Through  the  point  C  thus  obtained  draw  a  vertical  line  until 
it  intersects  line  AB.  Since  WC  is  the  horsepower  required  to  fly 
at  a  speed  equal  to  0.5  V  max.,  point  D  represents  a  plot  of  it 
drawn  at  the  proper  point  of  the  polar  curve.  Repeating  this 
construction  for  several  points  and  joining  all  these  points  by  a 
smooth  line  we  obtain  the  horsepower  required  plotted  against  the 
corresponding  lift  and,  by  virtue  of  construction,  also  the  known  speed. 


85 


Log  L 


FIGURE  5.  Determination  of  the  L/D  Ratio 


FIGURE  7,  Determination  of  Horsepower 
Required  at  Altitude 


Log  L 


FIGURE  6.  Determination  of  Horsepower 
Required  at  Sea  Level 


FIGURE  8.  Determination  of  Horsepower 
Available  at  Sea  Level 


HORSEPOWER  REQUIRED  AT  ALTITUDE 


To  obtain  the  curve  of  horsepower  required  at  any  altitude,  it 
is  sufficient  to  shift  the  sea  level  curve  to  the  right  by  the 
amount  equal  to  the  distance  (1)  which  is  measured  between  the  speed 
axis  and  the  ordinate  axis  at  the  height  corresponding  to  the  chosen 
altitude.  Figure  7  shows  this  construction  for  an  altitude  of  20,000 
ft. 


HORSEPOWER  AVAILABLE  AT  SEA  LEVEL 

To  construct  the  curves  for  horsepower  available  at  sea  level, 
we  introduce  three  new  constructional  curves  giving  the  variation  of 
thrust  horsepower  against  the  ratio  of  V/V  max. 

On  page  144,  Fig.  96,  of  Diehl’s  ’’Engineering  Aerodynamics,”  we 
fina  the  three  general  curves  giving  the  necessary  relationship.  In 
the  lower  left  part  of  the  diagram  we  incorporate  these  three  curves. 

The  scale  of  propeller  efficiencies  is  used  now  as  the  scale  of 
THP/THP  ,  while  the  scale  of  V/V  max.,  is  assumed  to  be  plotted  on 
the  axis  of  altitudes  (Fig.  8).  In  this  constitution  we  make  use  of 
the  line  (2)  from  Fig.  1.  The  power  plant  used  determines  whicn  of 
the  three  curves  shall  be  used.  The  construction  will  be  illustrated 
for  one  point  and  must  be  similarly  repeated  for  the  rest  of  the  points 
corresponding  to  fractions  of  V/V  max.  Figure  5  shows  the  construction 
for  point  of  V  equal  to  0.5  V  max.  See  also  Table  1. 

Having  chosen  one  of  the  THP/THPo  curves  we  project  the  corre¬ 
sponding  point  on  the  line  (2).  The  point  at  intersection  projected 
on  the  HP  scale  gives  the  horsepower  available  at  this  particular 
speed.  Drawing  line  (4)  unitl  it  intersects  the  horizontal  through 
the  point  on  the  polar  corresponding  to  the  speed  chosen,  we  obtain 
point  B  of  the  horsepower  available  curve  at  sea  level.  Repeating 
this  construction  for  several  more  points  and  joining  them  with  a 
smooth  curve,  we  obtain  the  curve  of  horsepower  available.  The 
curves  of  horsepower  available  and  horsepower  required  will  intersect 
at  point  B  of  Fig.  2,  corresponding  to  the  maximum  speed  when  the 
horsepower  available  is  equal  to  the  horsepower  required. 
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TABLE  1.  VARIATION  OF  THP  WITH  V/V  MAX  RATIO 


Ratio  of  THP /THP  Max. 

V/V 

max. 

R.P.M.  Max. 

1800 

2100 

2400 

1.2 

1.035 

1.015 

0.995 

1.1 

1.025 

1.015 

1.005 

1.0 

1.00 

1.00 

1.00 

0.9 

0.96 

0.97 

0.975 

0.8 

0.908 

0.925 

0.94 

0.7 

0.85 

0.865 

0.885 

0.6 

0.78 

0.80 

0.82 

0.5 

0.70 

0.72 

0.74 

0.4 

0.60 

0.617 

0.635 

0.3 

0.488 

0.495 

0.515 

HORSEPOWER  AVAILABLE  AT  ALTITUDE,  UNSUPERCHARGED 

To  obtain  the  horsepower  at  altitude,  we  will  incorporate  a 
curve  giving  the  variation  of  the  horsepower  with  altitude.  Since 
our  construction  will  be  made  assuming  that  the  speed  of  the  airplane 
remains  constant,  we  can  make  use  of  the  data  given  in  Table  10, 
page  139,  of  Diehl's  "Engineering  Aerodynamics".  See  also  Table  2. 

TABLE  2.  VARIATION  OF  THP  WITH  ALTITUDE 


Altitude 

THP/THP( 

Sea  level 

1.000 

5,000 

ft 

0.820 

10,000 

ft 

0.667 

15,000 

ft 

0.532 

20,000 

ft 

0.425 

30,000 

ft 

0.261 

32,000 

ft 

0.234 
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Figure  9  illustrates  the  construction  of  the  point  corresponding 
to  an  altitude  of  10,000  ft.  Inasmuch  as  the  transformation  of  the 
polar  for  any  altitude  is  made  by  an  Imaginary  shift  of  the  origin 
of  the  coordinates  along  a  45-deg  line  it  is  sufficient  to  locate  the 
point  B  on  the  horizontal  passing  through  the  10,000-ft  point  on 
the  axis  of  altitude  so  that  a  line  drawn  through  it  and  inclined 
45  deg  will  cross  the  axis  of  horsepower  giving  the  magnitude  CO 
equal  to  0.667  C A;  that  is,  if  the  point  A  corresponds  to  1000  HP, 
the  value  read  at  C  must  be  567  HP.  Having  the  curve  of  variation  of 
horsepower  with  altitude  it  is  sufficient  to  slide  the  curve  of  horse¬ 
power  available  at  sea  level  parallel  to  itself  by  the  amount  and 
along  the  direction  (1)  corresponding  to  the  desired  altitude  as  shown 
in  Fig.  9. 


FIGURE  9.  Determination  of  Horsepower  Available 
at  Altitude,  Unsupercharged 


HORSEPOWER  AVAILABLE  AT  ALTITUDE  WITH  SUPERCHARGED  ENGINE 

Most  of  the  present  supercharged  power  plant  installations  are 
made  so  that  the  original  sea  level  power  is  retained  up  to  a  certain 
fixed  altitude,  above  which  the  power  decreases  at  the  normal  rate. 

On  our  diagram  this  can  be  expressed  by  drawing  instead  of  the  previously 
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drawn  horsepower  variation  curve  a  new  curve  consisting  of  two 
branches:  a  vertical  straight  line  from  point  A  up  to  the  horizontal 
passing  through  the  altitude  up  to  which  the  sea  level  power  is 
preserved  and  a  curve  from  point  C  which  is  obtained  by  shifting 
horizontally  the  old  curve  until  points  B  and  C  coincide.  In  Fig.  10 
the  power  was  assumed  to  be  retained  up  to  an  altitude  of  10,000  ft. 
As  the  altitude  performance  on  the  polar  diagram  is  obtained  by  an 
imaginary  shift  of  the  origin  of  coordinates  along  a  45-deg  line, 
to  obtain  the  horsepower  available  at  10,000  ft,  we  must  shift  the 
curve  of  horsepower  available  at  sea  level  to  the  right  along  a 
45-deg  line  by  the  amount  equal  to  the  segment  (1). 


FIGURE  10.  Determination  of  Horsepower 
Available  at  Altitude,  Super¬ 
charged 


FIGURE  11.  Determination  of  Performance 
at  Altitude 


PERFORMANCE  AT  ALTITUDE 

To  obtain  the  speed  and  corresponding  horsepower  required  and 
available  at  altitude,  we  must  bear  in  mind  that  by  virtue  of 
construction  that  instead  of  a  horizontal  line  passing  through  W 
(as  it  was  done  in  Fig.  6)  it  is  necessary  to  use  another  horizontal 
line  such  that  if  W3  is  parallel  to  the  altitude  axis,  the  segment 
l-’B  must  be  equal  to  the  segment  OA  where  A  is  the  desired  altitude. 
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Figure  11  illustrates  this  for  an  angle  of  attack  corresponding  to 
the  maximum  speed  at  10,000  ft.  Since  the  point  of  intersection 
of  horsepower  required  and  available  curves  corresponds  to  the  max¬ 
imum  speed,  project  point  E  on  line  BC  and  OF.  Segment  OF  gives 
the  magnitudes  of  horsepower  required  and  available,  since  in  this 
case  they  are  equal.  Drawing  line  CD  parallel  to  speed  axis,  we 
obtain  the  magnitude  of  the  maximum  speed  given  by  the  segment  CD 
and  point  D  permits  us  to  interpolate  the  angle  of  attack  corresponding 
to  the  maximum  speed  at  10,000  ft. 

COMPLETE  PERFORMANCE  PREDICTION  WITH  THE  AID  OF  LOGARITHMIC  POLAR 
DIAGRAM 

From  the  foregoing  we  are  now  able  to  obtain  tne  complete  data 
on  theoretical  performance  of  a  new  design.  We  can  ascertain  the 
maximum  and  stalling  speeds  for  any  altitude  at  which  the  plane 
is  capable  of  flying.  By  joining  the  points  of  maximum  speeds 
plotted  on  the  corresponding  horizontal  lines  and  noting  the  inter¬ 
section  of  this  new  curve  with  rhe  line  tangent  to  the  polar  and 
parallel  to  the  speed  axis  (line  used  to  obtain  the  minimum  power 
required  for  flight),  we  can  determine  the  absolute  ceiling. 

To  obtain  the  best  speeds  for  climbing  and  to  determine  the 
excess  horsepower,  we  can  make  use  of  the  following  observation. 

From  the  analyses  of  a  number  of  airplanes  it  can  be  found  that  the 
best  speed  at  which  to  climb  at  sea  level  is  usually  equal  to  0.58 
V  max.  At  the  ceiling,  the  planes  would  fly  at  a  speed  corresponding 
to  the  mini, r. urn  power  necessary  to  maintain  flight.  If  we  join 
by  a  straight  line  the  point  corresponding  to  the  speed  equal  to 
0.58  V  max.,  on  the  horsepower-required  curve  at  sea  level  with 
the  point  of  minimum  power  necessary'  to  fly  at  the  ceiling,  the 
points  of  intersection  of  this  line  with  horse power- required  curves 
will  define,  within  the  accuracy  cl  practical  requirements,  the  best 
speeds  for  climbing.  Knowing  these  speeds,  we  can  obtain  the  excess 
horsepower  available  for  climb. 
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We  can  resort  now  to  standard  analytical  methods  and  obtain  the 
curve  of  rate  of  climb  and  time  to  climb  to  different  altitudes. 

Figure  12  shows  an  example  of  the  use  of  logarithmic  polar  as  it  is 
employed  in  routine  practice. 

CONCLUSIONS 

The  particular  advantages  afforded  by  the  method  described  in 
the  article,  in  addition  to  those  already  mentioned,  consist  of 
compactness  of  the  record  and  a  ready  means  of  estimating  slight 
changes  of  design.  In  many  instances  lack  of  time  does  not  allow 
an  elaborate  revision  of  a  performance  estimate  and  in  such  a  case 
the  logarithmic  polar  diagram  furnishes  an  easy  and  rapid  estimate  of 
the  effect  of  changes  in  weight,  parasite  resistance,  etc.  To 
the  trained  eye  the  graphical  picture  of  the  polar  gives  other  valuable 
information  such  as  the  general  efficiency  of  the  design,  the  approach 
to  inefficient  gross  flying  weight  condition,  and  such. 
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